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We study the dynamical spin response of doped two-leg Hubbard-like ladders in the framework 
of a low-energy effective field theory description given by the SO (6) Gross Neveu model. Using 
the integrability of the SO(6) Gross-Neveu model, we derive the low energy dynamical magnetic 
susceptibility. The susceptibility is characterized by an incommensurate coherent mode near (n,n) 
and by broad two excitation scattering continua at other fc-points. In our computation we are able 
to estimate the relative weights of these contributions. 

All calculations are performed using form-factor expansions which yield exact low energy results 
in the context of the SO(6) Gross-Neveu model. To employ this expansion, a number of hitherto 
undetermined form factors were computed. To do so, we developed a general approach for the 
computation of matrix elements of semi-local SO (6) Gross-Neveu operators. While our computation 
takes place in the context of SO (6) Gross-Neveu, we also consider the effects of perturbations away 
from an SO(6) symmetric model, showing that small perturbations at best quantitatively change 
the physics. 



PACS numbers: 



I. INTRODUCTION 



Ladder materials 1 - bridge the great divide between the 
well understood one dimensional (doped) Mott insulators 
and their perplexing two-dimensional analogs. The de- 
sire to understand their fascinating physical properties^ 
as well as the expectation that an understanding of lad- 
der materials would culminate in the theoretical un- 
derstanding of doped, two dimensional Mott insulators 
has led to intense, sustained theoretical interest in the 
problem. Analytical 9 ! 10 : 11 ! 12 : 13 ! 14 : 15 ! 16 ! 17 ! 18 ! 19 ! 43 as well 
as numericalSfliSiiS 2 * 2 ^ 2 ^^ techniques have been used to 
deduce the zero temperature phase diagrams of ladder 
models. The dynamical properties are less well charac- 
terized. In fact, in contrast to the half-filled caseSi, the 
accurate determination of dynamical response functions 
for two-leg ladder lattice models at finite doping remains 
a challenging open problem. One approach to calculat- 
ing dynamical correlations at finite doping is to utilize 
numerical methods such as quantum Monte-Carlo tech- 
niques or exact diagonalization. The single-particle spec- 
tral function and the dynamical spin and charge suscep- 
tibilities for both undoped and doped Hubbard ladders 
were computed by quantum Monte-Carlo methods in Rcf. 
i|22T) and compared to analytical random-phase approx- 
imation calculations for a two-leg Hubbard ladder. In 
Ref. |2f|) the dynamical structure factor was computed 
by exact diagonalization of a 2 x 16 t-J ladder model. 

A different approach is to utilize field theory methods 
to describe the low-energy degrees of freedom. In Ref. 
(|13() the structure of various dynamical response func- 
tions for the half-filled case was determined in the frame- 
work of a low-energy description in terms of the SO (8) 
Gross-Neveu model^ 7 .. Detailed calculations utilizing the 
integrability of the SO (8) Gross-Neveu model were car- 



ried out in Ref. I[43|) . In the doped case it has been 
argued in Refs. and l|TTl) that two-leg ladders can 
be described by the SO (6) Gross-Neveu model (£> U(l) 
Luttingcr liquid at low energies. The Luttinger liquid 
describes the gapless charge degrees of freedom on the 
ladder while the SO (6) Gross-Neveu model governs the 
behavior of the gapped degrees of freedom in the spin 
and orbital sectors. A prominent feature in the magnetic 
response due to bound state formation was analyzed in 
the framework of this theory and compared to numerical 
computations for a two- leg t-J ladder in Ref. |25|). 

It is important to go beyond the existing calculations 
in order to address questions raised by recent inelastic 
neutron scattering experiments on doped cuprate two-leg 
ladders^ 8 -. Furthermore, the dynamical susceptibilities of 
doped ladders can be used to develop a theory for systems 
of weakly coupled ladders, which in turn can be used as 
toy models for the striped phase in the cuprates^S. 

In the present work we determine the dynamical struc- 
ture factor of doped Hubbard-like two-leg ladders at low 
energies in the framework of the SO (6) Gross-Neveu de- 
scription. We go beyond Ref. l|25l) in that (1) we con- 
sider all momentum space regions with non-vanishing 
low-energy response; and (2) we consider multiparticle 
contributions to the dynamical structure factor. This 
allows us to evaluate the contribution to the structure 
factor of the bound state identified in Ref. |25|) relative 
to other contributions. 

Given that we employ the SO (6) Gross-Neveu model as 
a low energy description of the doped ladders it is worth- 
while dwelling briefly on the nature of this description. 
The SO (6) Gross-Neveu model represents the attractive 
fixed point of an RG flow of a doped ladder with weak 
generic interactions 1 — The flow has been established to 
two- loop order—. We point out that the flow to the SO (6) 
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Gross-Neveu does not assume that the bonding and anti- 
bonding bands of the Hubbard ladder have identical ve- 
locities. In the work of Ref. 10), it is shown that differing 
Fermi velocities will renormalize towards one another (an 
effect present only at two loops). 

As the RG is perturbative in nature, the bare cou- 
plings must not be particularly strong. How strong they 
may be before the RG becomes unreliable is indicated 
in Ref. (0). In this work, the related problem of half- 
filled ladders was studied where the RG flow takes the 
ladders to a sector where their low energy behavior is de- 
scribed by the SO (8) Gross-Neveu model. There it was 
found that U/t (U being representative of the strength of 
some generic short ranged interaction) may be as large 
as 6 before the RG flow is a poor indicator of the low 
energy physics. At larger values of U, it was found the 
charge and spin/orbital sectors decoupled. This decou- 
pling is already present in the U(l) Luttinger liquid ® 
SO (6) Gross-Neveu model and so the RG analysis for 
the doped ladders may well be descriptive for U/t > 6. 

Beyond a question of whether the bare couplings are 
sufficiently weak, it may be asked whether the RG flow 
accurately represents the physics in a wider sense. Where 
things may go wrong is illustrated by an appeal to the 
U(l) Thirring modcl^. Such a model has a Lagrangian 
of the form 

C = iyfl,^ + \g\\{ ]z f + \g±[(j x ) 2 + (jy) 2 ], (1-1) 

where jg = ^^^<j a ^ and ^,-0 are Dirac spinors. The 
one-loop RG equations for this model indicate that there 
is generic symmetry restoration to an SU(2) symmetric 
model. However the RG equations are misleading in the 
sector 7r — \g±] > —gn > \g±] > 0. In this sector the 
model maps onto the sine-Gordon model^ with interac- 
tion, cos(/3<i>), where (3 is given by 

/3 2 = 8tt-8/!; fi = cos _1 [cos(g||)/cos(.g_L)]. (1.2) 

(3 completely describes the physics of the model. And 
while the two coupling constants, gn and g±, flow under 
the RG, (3 does not, being an RG invariant. Thus the 
dynamically enhanced symmetry seemingly promised by 
the RG is here illusory. 

This case is however special. f3 as an RG invariant 
in the U(l) Thirring model is reflective of the presence 
of an SU(2) based quantum group symmetry. The RG 
flow cannot restore the SU(2) symmetry in this case as 
the symmetry has merely been deformed rather than 
broken*. There is however no counterpart of this ex- 
otic symmetry in the case of the SO(6) Gross-Neveu 
model and so this scenario of false RG-induccd symmetry 
restoration will not appear here. 

While the flow to an SO (6) symmetric model seems 
reasonably robust, we take a more cautious approach and 
understand the RG flow merely to move the model to a 
more symmetric state. Thus we should think of the low 
energy sector of the Hubbard ladders as an SO (6) Gross- 
Neveu model plus smallish perturbations. But because 



the SO(6) Gross-Neveu model is gapped, perturbation 
theory about it is well controlled: small perturbations 
lead only to small changes in the physics. We show in 
Section III that corrections induced by perturbations to 
SO (6) Gross-Neveu can indeed be straightforwardly han- 
dled and lead only to quantitative changes in the results. 

To extract the dynamic magnetic susceptibility from 
the integrability of the SO (6) Gross-Neveu model, we em- 
ploy a so called "form-factor" expansion of the spin-spin 
correlation function. This expansion yields exact results 
for response functions below certain energy thresholds. 
At the origin of this feature of the form-factor expansion 
is the notion of a well defined elementary excitation or 
particle in an integrablc model. Particles in an integrable 
model have infinite lifetimes - integrable models are thus 
in a sense superior versions of Fermi liquids. The in- 
teractions between particles in an integrable model exist 
only in restricted forms. The scattering of N particles 
can always be reduced to a series of two body scatter- 
ing amplitudes. This reduction occurs because of the 
presence of an infinite series of conservation laws. This 
simplified scattering yields that generic cigenfunctions of 
the SO (6) Gross-Neveu model can always be thought of 
as some distinct collection of the elementary excitations. 
In an non-integrable model this would not be possible 
because of particle production and decay. 

The particle content of SO (6) Gross-Neveu is relatively 
simple. The elementary excitations are found in two 
four dimensional representations, the kink (4) and the 
anti-kink (4). These excitations represent the elemen- 
tary fermionic excitations of the ladder (as discussed in 
greater detail in Section II). An interesting feature of 
SO (6) Gross-Neveu is that two kinks or two anti- kinks 
can form bound states. These bound states transform as 
the six dimensional vector representation (6) of SO (6). 
A portion of this six dimensional vector carries the quan- 
tum numbers of an SU(2) triplet. It is this bound state 
that provides a coherent response near (w, tt) in the doped 
ladders. 

With the spectrum of SO (6) Gross-Neveu in hand, 
we can readily write down the form factor expan- 
sion of any two-point correlation function (imagi- 
nary time ordered at T = 0) in the model: 

G°(x,r) = (O|T(0(x,T)0t(O,O))|O) 

oo 

=£5>|0(*,O)|n; S „> e -^x 

(n;s„|0t(O,O)|O>, (r>0),(1.3) 

where E Sn is the energy of the eigenstate, |n;s n ), and 
s n marks all possible quantum numbers carried by the 
state. In the second line we have inserted a resolution of 
the identity between the two fields of the correlator. This 
resolution is composed of eigenstates (i.e. multi-particle 
states) of the SO(6) Gross-Neveu Hamiltonian. We have 
thus reduced the computation of any correlation function 
to the calculation of a series of matrix elements represent- 
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ing transition amplitudes between the ground state and 
multi-particle states. 

The intcgrability of the SO (6) Gross-Neveu model al- 
lows the computation of these matrix elements in prin- 
ciple. Each matrix element satisfies a number of alge- 
braic constraints which can be solved. However as the 
particle number in the cigcnstatc \n : s n ) increases, solv- 

I 



where e = ± for fields, O, that are bosonic/fcrmionic. 
The delta functions present in the above equation guar- 
antee that only eigenstates with energy, to, contribute to 
the spectral function. As each particle has a fixed mass, 
higher particle states will not contribute if the sum of 
their masses is greater than the energy, uj, of interest. 

In this article, we compute only the one and lowest two 
particle form-factors of appropriate fields. This allows us 
to compute the spin response exactly to an energy of 2y/2 
times the spin gap in the region of k-space near (tx, it). In 
other regions of the Brillouin zone we are able to obtain 
exact results to energies twice the spin gap. However in 
these latter regions, the spectral intensity is considerable 
smaller than that found near (tt, 7r). 

While our results arc exact up to certain energies, it 
is likely that in certain cases (in particular, our compu- 
tation of the spin response in the region of (tt, 7t)) they 
remain extremely accurate beyond these thresholds. As 
a rule, form-factor expansions have been found to be 
extremely convergent with only the first terms yielding 
a significant contributioni2iiLi2iiS. As one example, in 
Ref. (Jilt) the spectral function for the spin-spin correla- 
tion function in the Ising model was computed. The two- 
particle eigenstates produced a contribution 1000 times 
greater than the four-particle eigenstates and a contribu- 
tion 10 s times greater than the six particle contribution. 
For practical purposes then the two particle form-factor 
in this case gave "exact" results for energies far in excess 
of the two particle threshold. 

The outline of this paper is as follows. In Section 
II, we review the connection between doped Hubbard- 
like ladders and the SO (6) Gross-Neveu model2iiii4. In 
particular we outline the renormalization flow that takes 
generically interacting doped ladders to the SO (6) Gross- 
Neveu model. We also delineate the connections between 
the excitations and operators in SO (6) Gross-Neveu and 
their corresponding ladder lattice counterparts. 

In Section III, we present our results for the dynamical 
magnetic susceptibility. We consider this quantity for a 
variety of dopings ranging from 1% to 25%. We also com- 



ing these constraints becomes increasingly cumbersome. 
Fortunately, if one is only interested in the behavior of 
the low energy spectral function (i.e. the imaginary por- 
tion of the corresponding retarded correlator), it is only 
necessary to compute matrix elements involving a small 
number of particles. The spectral function is given by 



(1.4) 



pare the results derived from our field theoretic treatment 
to an RPA analysis as well as to non-interacting ladders. 
We find there are significant differences. In particular 
the RPA cannot reproduce accurately the most promi- 
nent feature in the spin response, the coherent mode near 

(7T,7r). 

In Section IV, we present our results on form-factors 
in the SO (6) Gross-Neveu. The form- factors that are 
needed for the computation were not all known prior to 
this work. In part this came about because the spin re- 
sponse involve operators in SO (6) Gross-Neveu which are 
not standard objects of field theoretic interest. These op- 
erators are unusual in that they are semi-local while not 
being SO(6) currents. To derive the matrix of elements of 
these operators it was first necessary to establish the al- 
gebraic constraints satisfied by these form-factors. While 
it was not the primary intent of this work, we developed 
a general scheme by which the constraints for all semi- 
local operators in SO (6) Gross-Neveu (and by straightfor- 
ward extension, SO(2N) Gross-Neveu with N odd) may 
be written down. Portions of this developments are tech- 
nical. The basics of the derivation of all the form factors 
may be found in Section IV with certain details related 
to the semi- local operators relegated to Appendix B. For 
those interested in only the end product, a summary of 
all the relevant form-factors may be found in Section IV 
H. 



II. FROM DOPED LADDERS TO THE SO(6) 
GROSS-NEVEU MODEL 

Here we will briefly review the connection between 
doped Hubbard ladders and the SO(6) Gross-Neveu 
mode l 9 i 13 i 14 . We discuss the reduction of doped lad- 
ders to their field theoretic equivalent together with how 
to understand the excitation spectrum and fields of the 
50(6) Gross-Neveu model in terms of the original ladder 
model. 



1 f 

--ImG${x, -iw + S) = \ (0|O(x,0)|n; s n )(n; s n \O\0, 0)|0>*(w - E S J 

-e(O|0t(0, 0)|n; s„) (n; s n \0{x, 0)|0)<f(w + E S J 
I 



4 



A. Relation of doped D-Mott Hubbard Ladders to 
the SO(6) Gross-Neveu model 

We consider a weakly interacting ladder with generic 
repulsive interactions. It is well established^ that the sys- 
tem exhibits algebraically decaying d-wave pairing. We 
begin with non-interacting electrons hopping on a ladder: 

Hp = ~y] (ta\ a (x + l)a la (x) + ta\ a (x + l)a 2a (x) 

x.a 

+ t ± a{ a (x)a 2a (x)+h.c.). (2.1) 

Here the ai/a\ are the electron annihilation/creation op- 
erators for the electrons on rung I of the ladder, a; is a dis- 
crete coordinate along the ladder, and a =T,| describes 
electron spin, t and t± describe respectively hopping be- 
tween and along the ladder's rung. 

The first step in the map is to reexpress the a's of Hq 
in terms of bonding/anti-bonding variables: 

Cj a = —j={ai a + (-l) j a 2a ) 7 (2.2) 

(noting .7 = 1 = ab and j = 2 = b) . With this transforma- 
tion, the Hamiltonian can be diagonalized in momentum 
space in terms of two bands. As we are interested in the 
low energy behavior of the theory, the cj a 's are linearized 
about the Fermi surface, ItFj'- 

Cj a ^ Cf{j a e Fj -)- CLj a c Fj , (2-3) 

where L, R corresponding to the right and left moving 
modes about the Fermi surface. With this Hq becomes, 

H = dx ^2 v fj [c' Rjct id x c R j a - c' Lja id x c L j a ] . (2.4) 
J ]a 

Away from half-filling, the Vpj will be generically differ- 
ent. However for small to moderate dopings, the differ- 
ence of the two Fermi velocities will be small. We thus 
set VFat = vpi, in order to simplify matters. 

The next step is to consider adding interactions to 
the system. We consider adding all possible allowed 
interactions to this Hamiltonian. To organize this 
addition, we introduce various SU (2) scalar and vector 
currents 

l[j = CPia^aa'Cpja' 5 lfj = —Cpi T {e(T) TT iCpjr'; 
J Fj = C Ua C P3<? ; J S = 2 C PiT°VT'CPjT' • (2.5) 

Here P = R, L, a are the Pauli matrices and e CT(T < is the 
e-tensor. The crystal-momentum conserving interactions 
divide themselves into forward and backward scattering 
terms: 

Hb = ^2 tfj JijJij - bfj 3^ ■ Jfj ; 

i,j=l,2 

tt _ fP jR jL _ ra -wR tL 

"I Jij J ii J jj J ij J ii ' J JJ 



Here / and b are the forward and backward scatter- 
ing amplitudes. From hcrmiticity and parity, we have 
bn = b 2 \ and fu = hi- Furthermore, in the same spirit 
of taking the Fermi velocities of the two bands to be the 
same, bn = b 22 . Thus we obtain six independent ampli- 
tudes. As we are working away from half- filling in the 
chains, we do not consider Umklapp interactions. 

To facilitate the analysis of these interactions, we in- 
voke a change to variables. We begin by bosonizing the 
c's: 

c P]a = Kjae^" , P = +, — = R, L . (2.7) 

Here Klein factors satisfying 

{Kj a , Kip} = 2SijS al3 . (2.8) 

In terms of these four Bosc fields, four new Bose fields 
arc defined (effectively separating charge and spin): 

<t>Pi = -(0pit + + 0P2T + 0P2i) ; 

(j) P2 = -(</>pif — <j>pi[ + <t>P2] — 0P2|) ; 

0P3 = -('/'pit — 4>pii — 0P2| + 0P2i) ; 
p 

<t>Pi = -^(<f>Pl]- + <t>Pl[ - 0P2T - <t>P2i)- (2-9) 

Note that 4>P4 has a relative sign between the right and 
left movers. This sign effectively masks the underlying 
50(6) symmetry of the original Hamiltonian. The first 
and second bosons describe charge and spin fluctuations 
respectively. The latter two bosons, ^3 and </>4, are as- 
sociated with less transparent quantum numbers (as we 
discuss shortly). 

From these chiral bosons, one can define pairs of con- 
jugate bosons in the standard fashion 

<Pi = 4>Ri + <f>Li ; 

9 t = <t>m-<f>Li, (2-10) 

which obey the commutation relations, 

[ip(x),9(x')] = -iATre{x' - x) . (2.11) 

where Q(x — x') is the Heaviside step function. 

In terms of these variables the free part of the Hamil- 
tonian can be written 

#o = ^^{(^ Q ) 2 + (^ Q ) 2 }. (2.12) 

a 

The momentum conserving interactions can be written 

as 



(2.6) 



h b + h f = -±-^Y,A a {(d x e a ) 2 - (d x ip a ) 2 } 

-W{ 2 cos((9 4 ) cos(6> 2 ) - (612 ~ 46 i 2 ) cos(<^ 3 ) cos(6> 4 ) 
+ cos(0 2 )(26r i cos(0 3 ) + 2/f 2 cos((^ 3 )) 
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- cos(6> 4 )0i 2 + Ab p l2 ) cos(0 3 ), (2.13) 

where the coefficients A a are equal to A 1 / 4 = ±2(6^ ± 
/f 2 ) and A 2/3 = -(6f 1 ±_/ 1 CT 2 )/2. 

For generically repulsive scattering amplitudes, the 
various couplings / and b flow to fixed ratios under the 

p ri9.13.14 . 

b p 12 = ^i 2 = /i P 2 = -J&ii=5>0; 
fi2 = Ki = 0. 

6fa = J^i2 = 2/f 2 = -26? x = -i&fj = <? > 0; 
/fa = =0- (2.14) 

With these values, the interaction Hamiltonian dramati- 
cally simplifies to 

#int = H B + Hp 
4 

a=2 
4 

-4 9 ^ cos((9 a )cos((9 b ) . (2.15) 

a^b=2 

This Hamiltonian leads to gapped behavior for bosons 0, , 
i = 2,3,4. As we show below by a fcrmionization, the 
model governing the gapped behavior is the SO (6) Gross 
Neveu model. On the other hand we see that the inter- 
action Hamiltonian does not involve the boson governing 
total charge, 9\ — p+ . This boson remains gapless and 
so takes the form of a Luttinger liquid. If we were instead 
to allow Umklapp interactions (that is if the two chains 
were not doped away from half-filling) , this boson would 
be similarly gapped and the governing model would be 
SO (8) Gross Neveu. 

With the interaction Hamiltonian above, the theory 
describing the boson, 9 P+ , is a K = 1 Luttinger liquid. 
We will however consider deviations in K away from 1. 
Such deviations might occur either if the RG flow termi- 
nated before reaching the above set of fixed ratios or if 
differences in the Fermi velocity of the two bands were 
taken into account. Using the case of half-filled chains 
(equivalent to SO(8) Gross-Neveu) as a starting point, 
it was argued that generically K < 1 when doping was 
introduced^. 

We now recast this Hamiltonian explicitly in the form 
of the SO(6) Gross-Neveu model. We refcrmionize the 
bosons 4>p a , a = 1, • ■ ■ ,4 via 

*p n = K a e^ Pa , a= I,-- - ,3; 

^ P4 = P K4 e* 4 , (2.16) 

where the Klein factors are given by 

K\ = K 2 |, K2 = Kl|, K 3 =K U , K 4 = K 2 [ ■ (2.17) 



We then find the free Hamiltonian can be written as 
H = f dx J2 [*L ifl ***a - *L*^*i»] , (2-18) 

where the Fermi velocity, vp, has been set to 1, while the 
interaction Hamiltonian can be written as 

4 

H mt = -2g[J2(i*L*na-itf Ra *La)} 2 . (2.19) 

a=2 

This is precisely the interaction piece of the SO (6) Gross- 
Neveu model. It will sometimes prove convenient to re- 
cast the theory in terms of Majorana fermions, r\p a - In 
terms of the Dirac fermions, ^p a , they are given by 

*Pa = -^(V2a~2,P + «??2a-3,p) , (a = 2, 3, 4) . 

(2.20) 

In this basis, Hi nt can be recast as 

H mt = 2gG a p b G L b , (a^b= 1,...,6), (2.21) 

where G P b = irjp a ripb is one of the 15 SO (6) Gross-Neveu 
currents. 

B. Excitation spectrum of SO(6) Gross-Neveu 
model 

The SO(6) Gross Neveu has a total of 14 different 
types of excitations 3 ^. These excitations are organized 
into two four dimensional spinor representations and one 
six dimensional vector representation. The latter corre- 
sponds to the Majorana fermions introduced in the pre- 
vious section. We denote the fermionic creation opera- 
tors for the vector representation by A\, a = 1, ... ,6. 
The excitations corresponding to the spinor representa- 
tions are kinks and represent interpolations of the fields 
Ob, b = 1,2,3 between minima of the cosine potentials 
in Ean. l|2.15|l . There are eight different types of kinks, 
which we will denote by A^ a . Here a is a "multi-index" 
of the form a = (±1/2, ±1/2, ±1/2). The eight kinks 
divide themselves into two different four dimensional ir- 
reducible representations, something discussed in greater 
detail in Section IV 

SO (6) is a rank 3 algebra and the eigenvalues Ni, i — 
1, . . . ,3 of the three Cartan generators can be used as 
quantum numbers for the excitations. The combination 
of Majorana fermions 

m a -i±ir)2a, (2.22) 

generates excitations with quantum number N a = ±1, 
Nb = 0, b =/= a. The quantum numbers carried by the 
kinks A^ a are directly encoded in a. If a = (cti, a 2 , 0:3), 
oti = ±1/2, then A^ a creates a kink with quantum num- 
bers, Ni = oti. 

The quantum numbers, Ni, are related to the physi- 
cal quantum numbers of the system. In particular these 
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are the z-componcnt of spin. S z , the difference in z- 
component of spin between the two bands, 1S12, and 
the "relative band chirality", P12, defined as P\2 = 
Nri — Nli — Nr2 + Nl2, where Npj is the number elec- 
trons in band j with chirality P. In Ref. Il3l it was 
found: 



(JVi = 1,0,0) 


(S z 


= l,Sl2 


= 0,Pia 


= 0); 


(0,JV 2 = 1,0) <-> 


(S z 


= 0,5 12 


= l,Pl2 


= 0); 


(0,0,^3 = 1) ^ 


(S z 


= 0, 5i 2 


= 0, P 12 


= 2). (2.23) 



We note that these three quantum numbers are indepen- 
dent of the total charge, Q. Charge excitations all occur 
in the gapless sector governed by the boson 61 . 

We can see that the vector representation of excita- 
tions corresponds to states of two electrons (particle- 
holes pairs) in the original formulation. For example, the 
fermion rji Hr]2 carries spin, S z = ±1 and no charge and 
so is a magnon excitation. The spinor representations 
(the kinks) in turn are closely related to single particle 
excitations as their quantum numbers are combinations 
of Ni/2 (in fact they are single particle excitations mod- 
ulo their charge component). 



C. Identification of Physical Operators in terms of 
SO(6) Gross-Neveu Fields 

In this section we provide a dictionary between the 
SO (6) Gross-Neveu model and the original fermions of 
the Hubbard ladders. As discussed previously, the kinks 
are closely related to single particle excitations and so the 
original electron operators. There are 16 kinks in total 
(counting both left and right movers) and sixteen electron 
operators, the c's and cr's (four for each of the four Fermi 
points). In terms of the four bosons 4>RIL\ = 4>R/Lp+ an d 
4>R/Lii i — 2, 3, 4, the electron operators are given by c's 
are as follows: 



c fl6T ~ e 5*fli e 5(te-te-*R4). 

C Rbl ~ e i4>Rl e i(-4>R2+<t>R3-<pR4). 

CRab] ~ e 2*«i e 2(^ + 0H3+0fl4). 

CRabl ~ e ^Ri e l(-te-0 H 3+^R4) ; 



(even chirality) 



C Lbl ~ e i<pLl e i(-<f>L2+4>L3 + ll>L4). 

C Lb] ~ e l^l e i(^2-<Ai3+^4). 

CLabl ~ e i<t>Ll e H-<t>L2-<t>L3-<t>L4). 

CLab] ~ e i^l e iWx,2+^3-^4). 



The lattice fermions are all given by a charge piece, 
e l ^ R1 , multiplied by a field creating an SO(6) kink, 
e i l (±4>R2±4>Rz±4>Ri) _ The first four fermions involve even 
chirality kink fields while the latter involve odd chi- 
rality kinks. We denote the group representation of 
the even chirality kinks by 4. Electrons transforming 
as even chirality kinks are those above in addition to 

( c L6i' c Lbv c Labi> c Laht)' that is the Hermitean conjugate 
of the above odd kinks. The odd chirality kinks transform 
under the 4 representation. Again, electrons transform- 
ing as such are the above as well as the Hermitean con- 
jugate of the above even kinks, (c Rb1 - , c Rb[ , c Rab1 - , c Rabl ) 
The lack of symmetry between left and right movers re- 
flects the parity sign in the definition of <j)pi in Ean. Q2.9|l . 

To construct the spin-spin correlation function we need 
to consider the description of fermion bilincars in the lan- 
guage of the SO (6) Gross-Neveu model. Fermion bilin- 
ears can be of two types: one involving fermions asso- 
ciated with kinks of different chiralities and one involv- 
ing fermions associated with kinks of the same chirali- 
ties. In the first case we want to consider fermion bilin- 
ears transforming under 4®4 = 1©15. We will want 
to focus on the fifteen dimensional representation, 15, 
which is no more than the set of SO(6) currents. In 
the second case we want to consider bilinears behaving 
as 4 (g) 4 = 6 © 10 (and similarly for 4 ® 4). The six di- 
mensional representation, 6, is associated with particle- 
hole excitations carrying unit quantum numbers such as 
the magnon. The particle-hole excitations correspond- 
ing to the 6 are formed by taking symmetric combina- 
tions excluding diagonal terms of the even electrons in 
Ean. (|2~23|l : 



C Rb]J Labi +CRa b] c\ bi . 



(2.25) 



(odd chirality). 



(2.24) 



Such symmetric fermion combinations behave antisym- 
metrically under SO(6). As will be discussed, this corre- 
spondence between the 6 and symmetric combinations of 
fermion bilincars is most easily seen upon bosonization of 
the fermionic expressions. We will denote the fields creat- 
ing these two particle excitations as ip a , a = 1, . . . , 6. On 
the other hand, the ten dimensional representation, 10, 
arises from antisymmetric combinations of particle-hole 
excitations, 



C Rb1 C Labl ~ C Rab]C Lb[ , 



together with diagonal combinations of fermions 

CRb]c\ H . 

We will denote the fields creating these excitations by 
■0[ abc U where a, b, c — 1, ... 6 are anti-symmetrized. The 
odd representation for these fields will facilitate the com- 
putation of their associated matrix elements^. 
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III. SPIN-SPIN CORRELATION FUNCTIONS 
OF DOPED LADDERS 

In this section we utilize the form factors we calculate 
in Section IV to determine the low-energy behavior of the 
dynamical structure factor of the doped two-leg ladder. 
For readers uninterested in the details of the form factor 
computations, we have summarized the results in Section 
IV H. 



To evaluate these correlators we first express in terms 
of the bonding/anti-bonding fcrmions (a = 1,2) 



2S a = C ]b C db+c\ ab C lab - {-l) a (c\ b C lab +c\ ab C lb ) , 

2S a = c\b Cl ib + c\ ab C Ub - (-l) a (cjb C Tafc+4a6 C Tb) » 

(3.4) 



A. Preliminaries 

In order to compute the spin response, we need to re- 
late the lattice spin-spin correlation function to corre- 
lation functions in the SO (6) Gross-Neveu model. The 
dynamical magnetic susceptibilities of the doped ladder, 
q Xl q y ), is given in terms of two different spin corre- 
lators: one along a leg of the ladder 

Xn{x,t) = (Si(x I t).Si(0 I 0)) rot 

= (Safot) ■ S 2 (0,0)) rct = X22{x,t), (3.1) 

and one involving cross leg correlations 

Xi2(x,t) = (SiOM)-S 2 (0,0)) ret 

= <S 2 (.T,i) • S!(0,0)) rct EE X2l(.T,t). (3.2) 

Here () re t denotes the retarded correlation function. The 
dynamical susceptibility of the ladder is then given in 
terms of the susceptibilities along the legs and rungs as 

x( w >k) = Xn(v,kx) +X22(w, k x ) 

+2xi2(u,k x )cos(k y ). (3.3) 



The correlators we need can then be written (using SU(2) 
invariance) as 



/CO 
-OO 

x(TS a (x,r)-S 6 (0)) 
= (S a -S b )(u,k). (3.5) 



The dynamical spin-spin correlation functions are ex- 
pressed as 

^Sx-SiXw.fc) = (S+S^)(u;,k) = ^(I 1 +I 2 + I 3 ); 

H(S 1 -S 2 )(w,fc) = (5 1 + 5 2 -)(c,fc) = i(/ 1 + / 2 -/ 3 ), 

(3.6) 

where the It are given by 



h(u;,k) = I drdxe^+^-^iTcl^rMx, tJc^ (0)^(0)), 
h{u,k) = J drfa i (- iu + 6 ) T - fe (T4 t ( a; ,T)c oH (x i T)4 H (0)c a6t (0)), 

h(v,k) = J dTdxe^ + ^- lk ^T[ci b ^ Xl T)c H (x,T)+^ [^(0)^(0) + 4 H (0)c feT (0)]). (3.7) 

In the next step we take the low-energy limit, in which the fermion operators are expanded according to l|2.3|l . This 
gives, for example 

4t ( x > T ) c bi ( x > T ) ~ c Rb] ( x > T )°Rbl (x,t)+ <4 bT (x, t)c lh (x, t) 

+e~ 2iKFbX J m (x, r)c LbL (x,t)+ e %iK ^ x c\^ (x, r)c Rbl (x, r) . (3.8) 

Substituting (|3.8(l and its analogs into the expressions for Ij(u>,k) we obtain sums of terms, each of which can be 
calculated using a form factor expansion as sketched in the introduction. For example in the imaginary piece of 
I\{ui, k) we have a contribution of the form 

OO ft Tt 

EE / II 1 )n d 2 d " I (oi4, T (° ■ ■ ■ 0nK...*„\* 5 ( u - E m ^ coshg ^-E?^). ( 3 - 9 ) 

n=0{ arl } 71 '' ' i=l i=l VF 



8 



where the first sum J2 n runs over the number of particles 
in a particular term of the form factor expansion and the 
second sum runs over the different possible particle 
types. The complete set of intermediate states used in 
the spectral representation (|3.8|) is defined by 

\e 1 ...e n ) ai ... an =Ai n (6 n )---Al(6 1 )\0) . (3.10) 

Carrying out the expansion 12.3|l for all terms in (|3.7() we 
obtain 



h{u),k) 



AhJi(lli, k) if k 

A 12 J 2 {uj,k^2K Fb ) if k 





±2K Fb , 

(3.11) 



A 21 Ji(u,k) iffcwO 
A 22 J 2 {uJ,kT'2K Fab ) \ikK±2K Fab , 

(3.12) 

j A 31 J!(uj,kT K F ^) iik^±K F - 
y A 32 J 3 (lu, k T Kf+) if k « ±K F+ . 

(3.13) 



Here we have defined 



K F ± = K 



Fab 



±K 



Fb ■ 



(3.14) 



The Ji's correspond to the lowest (one and two) parti- 
cle (bound state and kink) contributions to the spin re- 
sponse. Their explicit forms are derived in Appendix lAl 
At low energies, as explained in the introduction, these 
terms provide the sole contributions. The various A' s in 
the above are normalization constants. To fix them, we 
compare the spectral weight in the spin response obtained 
using our low energy field theoretic reduction to that 
obtained using a random phase approximation (RPA). 
While the RPA is a crude approximation, we believe that 
it will provide a rough estimate for the integrated spin 
response characterized by 



Imx(w, k) 



IV 



(3.15) 



Thus by equating the weight given by the two approaches, 
we are able to fix the values of the A's. The details may 
be found in Appendix lAl 



B. Spin Response at U = t 

In this section we discuss the spin response at the rel- 
atively moderate value of U = t. At the top of Figure 
1, we present an intensity plot from the field theoretic 
computation for the spin response function at 1% dop- 
ing, both at k y — (the left panel) and k y = it (the 
right panel). The low energy spectral weight is found 
near values of k x related to K Fab and K Fbl the bonding 



and anti-bonding wavevectors. As a function of doping, 
jj,, away from half-filling, the values of K Fab and K Fb are 



K Fb = cos 



' 2t 

.2 



r 2 



t 2 j_/t 2 - sin 2 (7172), 1/2' 

1 + COs(7T/i) ' 



K Fab = 7r(l - n) - K Fb . 



(3.16) 



At k y = 7r, the low energy weight is present near k x — 
K Fa b + K Fb , K Fb - K Fabl 2ir - K Fab - K Fb , and 2n - 
K Fb +K Fab . We see that the greatest intensity is found at 
k x = K Fa b + K Fb and k x = 2tt — K Fb — K Fab - Exactly at 
half-filling, these wave-vectors equal 7r, i.e. K Fab +K Fb = 
it. Thus as the system is doped the (it, it) peaks splits 
into two incommensurate peaks. However at 1% doping, 
the splitting is barely resolved in Figure 1. For k y = n, 
the spectral weight at k x = K Fb — K Fab is considerably 
smaller. 

Spectral weight is also found for k y = at k x = 2K Fabl 
2ir — 2K Fab , 2K Fb , 2ir — 2K Fb , and finally at k x = 0, 2ir. 
As we can see from the intensity scales, it is considerably 
less than that at k v — ir. We also see that the energy at 
which weight is first found is higher. This reflects the fact 
that the first contribution to the spectral weight is a two- 
particle contribution and so is at an energy a y/2 times 
higher than that at k y = ir where the initial contribution 
is made by a single coherent bound state. 

In the lower and middle panels of Figure 1 , we compare 
our field theoretic results for the spin response function 
both to the results for non-interacting ladders (see Ap- 
pendix C) and to interacting ladders treated using an 
RPA approximation (see Appendix A 4). The RPA ap- 
proximation sees a more diffuse distribution of spectral 
weight than that of the field theoretic treatment. In par- 
ticular the RPA approximation does a poor job at deduc- 
ing the existence of the contribution of a single bound 
state to the spin response near (ir, ir). But both the RPA 
and field theoretic treatment predict a much larger low- 
energy response than that of non-interacting ladders. In 
this sense the ladders mimic the behavior of individual 
chains where the presence of interactions shifts spectral 
weight from high to low energieat The presence of low 
energy spectral weight in a single chain is captured by 
the Miiller ansata&. It would be interesting if a similar 
ansatz might be developed for the spin response of doped 
ladders using our field theoretic treatment as a starting 
point. 

To obtain a more quantitative understanding of the 
spin response, we have plotted constant wavevector cuts 
of the spectral function in Figure 2. Here we have done 
so at 10% doping so that the low energy spectral weight 
present near any of the values of k x related to K Fab or 
K Fb docs not overlap with weight arising from a dif- 
ferent such value of k x . We first consider the cut at 
K Fa b + K Fb . We see explicitly the onset of weight at 
the spin gap, A s = \/2m, the excitation energy of the 
bound state. We note that contribution of the bound 
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Non-interacting 



FIG. 1: Plots of the spin response function computed in three different fashions: using SO(6) Gross-Neveu as a low-energy 
effective theory; RPA; and taking the ladders as non-interacting. The left-hand panels give the response at k y — while the 
right-hand panels correspond to k y — n. In making these plots, we have assumed 1% doping and U — t. We have plotted k x 
in rationalized units of 2ir. 
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FIG. 2: A number of constant wavevector cuts of the spin response function for 10% doping. We compare the computations 
done in the SO(6) Gross-Neveu framework with both non-interacting ladders and an RPA treatment. 
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state does not lead to a (theoretically) infinitely narrow 
feature in the spin response. While this would be so at 
half-filling, the gapless nature of the charge mode at finite 
doping smooths the feature out giving it a finite spread. 
Nonetheless there is a divergence at (to 2 — k 2 / vp) — > go- 
ing as (uj 2 — k 2 /v 2 F ) K / 2 ^ 1 . In the energy interval (0, 10m) 
we can also estimate the amount of spectral weight due to 
the single bound state. We find that it accounts for 37% 
of the total spectral weight (with a Luttinger parameter, 
K = 0.93). Of this 37%, 42% of it occurs in the interval 
(0,2m), i.e. before the incoherent two kink contribution 
to the spectral function begins to be felt. 



C. Spin Response at Larger U 

As U becomes larger we expect that the Fermi velocity 
will renormalize significantly but that our effective low 
energy description will remain valid. This view is sup- 
ported by the two loop RG computations of Fabrizio^. 
There it was shown that the two-loop RG preserved the 
flow to SO (6) but that unlike one- loop RG equations, the 
Fermi velocity was rcnormalizcd downward. 

To study the spin response at large U we thus suppose 
that the ratio Vf/tti is much smaller than at weak cou- 
pling. The RPA breaks down at U in excess of t, and so 
we cannot directly estimate the corresponding normal- 
ization constants, A's. Instead we simply continue to use 
the values calculated for U = t. We point out that while 
at best any comparison of low energy spectral weight at 
different wavevectors will then be qualitative, the spec- 
tral weight at any one value of k will be unaffected by 
this procedure. 

In Figure 3 we plot the spin response for a value of U 
corresponding to vp/m — 4 and for three different values 
of doping, 5%, 10%, and 25%. In comparison with the 
spin response in Figure 1 where vf/tu ^> 1, the response 
at a given energy occurs over a much wider band of k x 
values. 

From Figure 3, we can observe how the spin response 
evolves as a function of doping. Two primary effects are 
observed. The first is increasingly significant incommcn- 
suration effects as doping is increased. Incommensura- 
tion effects are most noticeable near {%, tt), the region 
of greatest spin response. At small doping the sum of 
bonding and anti-bonding wavevectors remain close to tt 
(see Ean. H3.16jl '). Consequently at small dopings the spin 
response near (71", tt) appears akin to that at half-filling 
with only a single peak of intensity apparent. As dop- 
ing is increased, however, two regions of intensity, one at 
KFab+Kpb and one at iTT—Kpah—Kpb, become resolved. 
At 5% doping, the two regions are becoming individually 
visible. As we increase doping further, the regions sep- 
arate and become distinct over a wider and wider range 
of energy. In Figure 4, we plot the spin response in a 
narrow region about (tt, tt) for the three different values 
of dopings. In this figure, the splitting of the (tt, tt) peak 
as function of doping is more readily apparent. 



An obvious signature of the incommensuration in a 
neutron scattering experiment would be found in a series 
of constant k-scans as a function of energy. In comparing 
such scans at (tt, tt) and (Kp a b + Kpb,Tr), a greater in- 
tensity at the latter value of k x would mark the presence 
of incommensuration. 

Incommensuration effects will be less obvious at k y = 
simply because the intensity present at this value of 
k y = is markedly less than at k y = tt. But the k x 
points at which low energy intensity is found, i.e. k x = 
2Kp a b, 2Kpt, and will shift relative to one another. 
This effect is convoluted with changing relative intensity 
at k x = 2Kpb and k x = 2Kpab- For smaller dopings 
(i.e. 5% and 10%), the intensities found near these two 
points is roughly equal. For larger dopings, however, the 
intensity at 2Kp a b becomes notably larger. In the Figure 
3 panel for k y = at 25% doping, the contribution at 
k x = 2Kp a b and k x = overlap, while the intensity at 
k x = 2Kpb is not observable in the plot at hand. 

The second notable effect of doping is that with in- 
creasing doping one sees a decrease in maximum inten- 
sity at ky = tt relative to that at k y = 0. At 5% doping, 
the relative intensity at k y — tt is twelve times that at 
k y = 0. At 25% doping, this ratio has been reduced from 
twelve to two. 



D. Effects of Integrability Breaking Perturbations 

In this section we consider the effects that integrabil- 
ity breaking perturbations will have upon our results. 
As we discussed in the introduction, this is important to 
consider because the renormalization group flow leaves a 
generic system more symmetric but does not necessarily 
render the low energy sector exactly equivalent to the 
SO (6) Gross- Neveu model. As was done for the SO (8) 
Gross-Neveu description of half-filled ladders^, we con- 
sider this question in broad terms so as to demonstrate 
that small perturbations away from the SO (6) Gross- 
Neveu model do not drastically change the results. (For a 
more detailed analysis for general ladder models see Ref. 

©•) 

With this in mind we examine general possible per- 
turbations of a U(l) Luttinger liquid (g> SO(6) Gross- 
Neveu model. Possible perturbations can be read off 
from Ean. H2.13|) . the Hamiltonian of the doped ladders 
with generic (bare) interactions. These interactions all 
take the form of 

# P crt = Acos(6>i) cos(^), % 7^ j (3.17) 

or 

H peit = Xcos((p 3 )cos(6j). (3.18) 

In general these interaction preserve the charge, spin, and 
relative band chirality but not the difference in the z- 
componcnt of spin between bands (owing to the presence 
of cos(<£ 3 )). 
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FIG. 3: Plots of the spin response function for three different dopings, 5%, 10%, and 25%. The left-hand panels give the 
response at k y — while the right-hand panels correspond to k y — n. These plots are made assuming a strongly renormalized 
Fermi velocity, vf/uti = 4. 
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To determine the manner in which the above type of 
perturbations will affect the spectrum of the model is 
straightforward. We can immediately note that the var- 
ious multiplets of SO (6) Gross Neveu will be split by 
the perturbation, but the single particle multiplets them- 
selves will not be mixed. We expect also that the notion 
of single particle excitations will remain robust - no per- 
turbation in Ean. (|3.16[) will lead to kink confinement. 
Moreover no perturbation will gap out the U(l) Luttingcr 
liquid. 

To determine the splitting of any given multiplet of 
particles {Ai(0)}, we employ stationary state perturba- 
tion theory. In this sense we operate in the same spirit 
as Ref. (|31|) in treating the off-critical Ising model in a 



magnetic field or Ref. (|32f) in treating perturbations of 
the 0(3) non- linear sigma model. At leading order this 
amounts then to diagonalizing the matrix 



Mi. 



(M8)H pcit A}(0)) 



(3.19) 



The computation of these matrix elements is straightfor- 
ward in the context of integrability. Each matrix element 
is no more than a form factor which can easily be com- 
puted in the same fashion as done in Section IV. More- 
over because the theory is massive, perturbation theory 
is controlled and so small perturbations should only in- 
troduce small quantitative changes to our results. 

Deviations in the mass spectrum will have a two fold 
effect. The bound state producing the peak in the spin 
response at uo = \/2m and k = KFab + KFb will, un- 
der the perturbation, be mixed into other members of 
the vector multiplet, leading to a number of particles 
that will couple to the spin operator. Consequently the 
peak will split. This splitting however will be obscured 
by the gapless charge excitations with their tendency to 
smear any sharp spectral features. This smearing will be 
enhanced by finite temperature. Unless the splitting is 
then large, it is unlikely the effects of the perturbation 
will be detectable experimentally 

Like the splitting of the multiplet of bound states, the 
two particle threshold will similarly split. Rather than 
occur at u> = 2m, a number of thresholds will appear 
about this energy. This effect will be muted at points in 
k-space where the threshold opens smoothly from zero 
rather than appearing as a singularity. But the effect 
will be more pronounced at k x = 0, k y = and k x = 
Kpb — Kpab, k y = 7r, where the two particle threshold is 
marked by a van-Hove like square root singularity. How- 
ever the relative weight of spectral intensity is small at 
these points (in comparison to k x = Kpab + Kph: k y = tt) 
and so even here any splitting will presumably be difficult 
to detect. 

We can also consider the effects of integrable break- 
ing perturbations upon matrix elements. The behavior 
of matrix elements near threshold determine whether a 
threshold opens up smoothly (for example the two par- 
ticle matrix element vanishing at threshold, i.e. {uj 2 — 



k 2 /vp) — > 0), extinguishing a van-Hove singularity) or 
whether it is marked by a discontinuity (the matrix el- 
ement is finite at threshold). Small perturbations may 
then lead a vanishing matrix element to be finite forcing 
a qualitative change in the physics. However, we now 
argue that this does not happen at leading order in the 
perturbation. 

To compute the leading order correction of a matrix 
element of some operator, 0(x = 0, t = 0), we employ the 
relation (arrived at from Keldysh perturbation theory) 



(S(0|O(0,0)|i) = 

dt / dx(Q\[Q{Q),H peTt {x,t)]\i)- (3-20) 

-oo J — oo 



Here |0) is the vacuum state and \i) is some other state 
with some arbitrary number of particles. If we employ a 
form-factor expansion to evaluate the expectation values 
involving operator bilinears in the above expression we 
obtain 



(5(0|C>(0,0) |i 

M coo 



dt f dxY,((0\O(0)\j)(j\H pci . t (x,t)\t) 

oo J — oc j 

(0\H pelt {x,t)\j)(j\O(0)\i)), (3.21) 



where ■ represents schematically a sum over all states 
appearing in a resolution of the identity. We argue be- 
low that this correction will uniformly vanish for ma- 
trix elements, involving two kinks or anti- kinks with 
minimal energy, 2m, and for which the matrix element, 
(0|O(0, 0)|i), itself vanishes. This implies that the lead- 
ing correction leaves unaffected the behavior of all rele- 
vant matrix elements at threshold appearing in the spin 
response. 

In the relevant cases where the matrix element itself 
vanishes, the vanishing is a consequence of the structure 
of the low energy S-matrix. The S-matrix at vanishing 
energies of two kinks is — P, where P is the permuta- 
tion matrix. Not only then docs (0|O(0, 0)|i) vanish, 
but any matrix element involving \i) vanishes, for exam- 
ple (j\H pert (x,t)\i) and (0\H pcit (x,t)\i). The vanishing 
comes about through an application of the scattering ax- 
iom (see Section IV B). This axiom implies the following 
relation, 



(j\H peit (x,t)\i) = -(j\H peIt (x,t)\i). (3.22) 

This then necessarily implies the matrix element vanishes 
regardless of the nature of \j) or H pelt . 
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FIG. 4: Plots of the spin response function for three different dopings, 5%, 10%, and 25% for k y — n and a narrow range of k x 
chosen to emphasize the region where the maximum intensity is found. These plots are made assuming a strongly renormalized 
Fermi velocity, vp/m — 4. 



IV. FORM FACTORS FOR THE SO (6) 
GROSS-NEVEU MODEL 

A. S-matrices 

In this section we lay out the S-matrices for scattering 
between the kinks. From these, we will be able to develop 
formulae for the various form-factors needed to compute 
the spin-spin ladder correlation functions. 



1. SO(6) Spinor Representations 

The spinor representations are expressed in terms of 
SO(6) 7-matrices. The 7-matrices in turn are built out 
of the two-dimensional Pauli matrices, <r l 's. As we are 
interested in 5*0(6), we consider three copies of the <7 l 's, 

4,4,4, (4-i) 

each acting on a different two-dimensional space, i.e.: 

a l 2 \ai,a 2 ,a 3 ) = \ai,ala 2 ,a 3 ). (4.2) 

A particular basis vector, \a), in the corresponding 8 di- 
mensional vector space can be labeled by a series of three 
numbers ±1/2, i.e. 



\a) = K,a 2 ,a 3 ) = | ± 1/2, ±1/2, ±1/2). 



(4.3) 



Physically, the kink associated with \a) then carries 
quantum numbers Ni = ±1/2 for the three U(l) sym- 
metries corresponding to the Cartan elements of SO (6), 
discussed in the previous section. 

In terms of the a' 1 's, the 7-matrices are defined by (fol- 
lowing the conventions of Ref. l33|) , 



72n-l 



n 

*!=i 



l2r. 



n 

*!=i 



(4.4) 



for 1 < n < 3. These matrices satisfy the necessary 
Clifford algebra: 



{7a, lb} = 2S ab . 



(4.5) 



In this representation the Clifford-algebra generators, j n , 
are imaginary and antisymmetric for n even, while for n 
odd, they are real and symmetric. The 50(6)-gencrators 
are represented by 



7V), 



(4.6) 



in analogy with the more familiar 50(4) case. 

The 8-dimensional space of the 7's decomposes into 
two 4-dimcnsional spaces, each of which forms one of the 
two irreducible 50(6) spinor representations. The de- 
composition is achieved explicitly via the hermitian chi- 
rality operator, T, 



r = H) 3 II 7 a = °*i 



(4.7) 



r is such that it commutes with all 50(6) generators and 
is diagonal with eigenvalues ±1. If | ct ± ) is a state with an 
even (odd) number of negative components (i.e. states 
with positive or negative isotopic chirality), then 



(4.8) 



Thus the operators (1 ± T)/2 project onto the two irre- 
ducible subspaces. When it is necessary to make a dis- 
tinction, we will denote kinks with an even number of 
negative components by a and those with an odd num- 
ber by a. 

The last item to be presented in this section is the 
charge conjugation matrix, C . In terms of the 7's, C is 
given by 



C = 717375- 



(4.9) 
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C is completely off-diagonal (as expected). Moreover we 
note that C is skew symmetric (C T = — C, C 2 = —1). 
This fact will have important consequences in the deter- 
mination of the form factors. 

If a = \ai,ct2,a3) is an even kink its corresponding 
anti-particle is denoted by a = | — ai, — ot2, — ots). By the 
skew-symmetry of C, we thus have 



\a) = (-ir 2 - 1/2 C a c 



2. S-Matrices for the Kinks 



(4.10) 



In order to describe factorizable scattering, we intro- 
duce Faddeev - Zamolodchikov (FZ) operators, A^ a (6), 
that create the elementary kinks. 9 is the rapidity which 
parameterizes a particle's energy and momentum: 



, r(l/4-0/2wr) 
T(l/4 + (9/2i7r)' 



(4.16) 



P{9) =msinh(0); 



E{9) = mcosh(6>). (4.11) 



We note that there is a sign difference in the definition of 
u\ between Ref. I33I and Ean. I|4.16|) . This sign difference 
arises as we give here the S-matrices of the physical par- 
ticles (those with fractional statistics). S^(9) satisfies a 
Yang-Baxter equation, 

si\l\ (o 12 )s^n (M^ = 

(M^ (0i2)- (4.17) 

Physically, the Yang-Baxter equation encodes the equiv- 
alence of different ways of representing three-body in- 
teractions in terms of two-body amplitudes. Formally, it 
expresses the associativity of the Faddeev-Zamolodchikov 
algebra. The S-matrix, S^g{S), also satisfies both a cross- 
ing, 



Because the 50(6) Gross-Neveu model is integrable, 
scattering is completely encoded in the two-body S- 
matrix. This S-matrix, in turn, is parameterized by the 
commutation relations of the Faddeev-Zamolodchikov 
operators: 

4(0i)4(0 2 ) = Sl${9 12 )A\{6 2 )A\{0i)- (4-12) 

^12(^12) is the amplitude of a process by which particles 1 
and 2 scatter into 3 and 4. It is a function of 9 12 = 9\ — 92 
by reason of Lorentz invariance. 

Kink-kink scattering takes the form 



r=0 



°7/3 a Sa ' 



(4.13) 



where is a rank-r antisymmetric tensor: 
a y ' = a =(7 7 )a- 
Here A represents a complete anti-symmctrization of the 

(r) (r) 

gamma matrices. By cr^cr^ we mean a trace over all 
possible rank-r antisymmetric tensors: 



(4-14) 



_M (r) 
a <yl3 a Sa 



E *?f ar *r a " ar - (4.i5) 



The generic form of S^fa was determined in Ref . l36l while 
the specific forms of the it's were given in Ref. [33. There 
it was found 



u 3+r (9) 
u r+2 {9) 
u 2 {9) 



(-l) r u 3 - r (9); 

, a , (l-r/2)-(l + i0/n) 
M >(l-r/2) + (l + i9/ny 
T(l - 9/2iir)T(8/2iir -1/4) 
r(3/4-6»/2i7r)r(6»/2i7r) : 
3/2- 9/inr(l/2 + 9/2iTr) 
1 -9/in r(l/2-6»/2i7r) 



S^{m-9) = S^{9)C^C a& . 
and unitarity relation, 



-,7(5 
'a/3 



(o)s;'f'(-8) 



"a °3 



(4.18) 



(4.19) 



We note that for the above kink-kink S-matrices, crossing 
is satisfied, as expected, with C = 717375. 

These constraints determine SZ. S g up to a 'CDD'-factor. 
Such factors allow additional poles to be added, if needed, 
in the physical strip, Re(6*) = 0, < Im6> < 2ir, to the 
scattering matrix and so are indicative of bound states. 
Here we expect to find a pole in u n , n even, at 9 = 
in/2 reflecting the fact that two (same chirality) kinks 
can form a bound state (the fundamental fcrmions of the 
model) of mass \/2m. The overall sign of the S-matrix is 
determined by examining the residue of the pole in S^g 
at 9 = iir/2. This pole is indicative of the formation of 
a mass y/2m bound state in the s-channel and so should 
have positive imaginary residue. 

The 50(6) Gross-Neveu model has an isotopic chirality 
conservation lawS. Thus opposite chirality kink scatter- 
ing is determined solely by u n , n odd, while same chirality 
kink scattering is determined solely by u n , n even. For 
example, it can be shown that for even chirality kinks, 
the S-matrix reduces to 



S 2p(0) = -u 2 {9)5 ia 5 S B 



u 2 {9))5 ia 8 Sa . (4.20) 



B. Basic Properties: Two Particle Form Factors 

The two particle form factors of a field ip(x) arc defined 
as the matrix elements, 

ff 2 (9 1 ,9 2 ) = (OWO)4(0 a )4(0i)|O). (4.21) 
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The form of f 12(61:62) is constrained by integrability, 
braiding relations, Lorentz invariance, and hermiticity. 

The constraint coming from integrability arises from 
the scattering of Faddeev-Zamolodchikov operators. For 
Ean. (|4.21|> to be consistent with Ean. l|4.12[> . we must 
have 

ft 1 (62:6 1 ) = Slt(6 12 )fl(6i:62). (4.22) 

The second constraint can be thought of as a periodicity 
axiom. It reads 

fti(62:6i) = Qi2R^ li) ff 2 (6i - 2™, 62). (4.23) 

Both Q12 and are phases. Q12 arises from the 

charge conjugation matrix in SO (6) Gross-Neveu being 
skew-symmetrio^l. Its determination is discussed in some 
detail in Appendix B. The phase factor R^^, on the 
other hand, arises from the "semi-locality" of the fields 
in the 50(6) Gross-Neveu model. Working in Euclidean 
space, fields can be said to be semi-local if their prod- 
uct sees a phase change when the fields are taken around 
one another in the plane via the analytic continuation 
z — > ze l2jr , z — > ze~ t27T . More specifically we have 

V>i(ze 42T , ze~ i2n )ijj(0) = R^M*: z)V>(0). (4.24) 

Here ipi is the field that is associated with the particle 
A\(Oi). The phase R^^p can be calculated from the op- 
erator product expansions of the operator ip with bosonic 
vertex operator representations of the kink fields. If A\ 
is a kink, \a) = (±1/2, ±1/2, ±1/2), the corresponding 
representation of the field is 

fa = e H±0i+±<fe+<fe) ( 4 25) 

In contrast to the representation of the left and right 
components of the lattice fermions fEan. H2.24fl ). these 
fields are non-chiral. The form factor, ff 2 , must also 
satisfy Lorentz covariance. If ip has Lorentz spin, s, ff 2 
will take the form (at least in the cases at hand), 

ff2(6i:62) = e s ^ +e ^ 2 f(6 12 ), (4.26) 

where by virtue of 612 = 61—82, / is a Lorentz scalar. 

The constraints Ean. (|4.22[> . Ean. (|4.23|l . and 
Ean. (|4.26|l . do not uniquely specify ff 2 - It is eas- 
ily seen that if ff 2 satisfies these axioms then so docs 
fi2 fl , ( c °sh(^i2)), where g(x) is some rational function. 
The strategy is then as follows. One first determines the 
minimal solution to the constraints, minimal in the sense 
that is has the minimum number of zeroes and poles in 
the physical strip, Rc(8) = 0,0 < Im(6>) < 2tt. Then 
one adds poles according to the theory's bound state 
structure. If the S-matrix element describing scattering 
of particles 1,2 has a pole at 8 = iu, then ff 2 (6i, 6 2 ) has 
a pole at 

8 2 =6 l +iu. (4.27) 



Insisting that f{ 2 has such poles and only such poles fixes 
g(x) up to a constant. 

The phase of this constant can be readily determined. 
Appealing to hermiticity gives us 

(0^(0)4(^)4(^)10)* = (01^(^)^(02)^(0)10) 

= C I1 C2 2 (O|V t (O)4(0i - ^7^)4(0 2 - wr) 1 0) ,(4.28) 
where the last line follows by crossing and so 

ff 2 (8i,6 2 y = 0-^0-22^(62 - in, 8 X - in). (4.29) 

For the kinks the matrix C equals the standard charge 
conjugation matrix, C = 717375, while for the Majo- 
rana fermions (7 = 1. (We note however that C for the 
Majorana fermions is not identically 1. This difference 
between C and C for the Majorana fermions arises be- 
cause under Hcrmitcan conjugation these fermions are 
invariant while under charge conjugation they undergo 
sign changes.) For the form factors wc will examine, her- 
miticity will be enough to fix their overall phase. 

As a final consideration we examine the large rapidity 
(8 — > 00) asymptotics of a form factor. These asymp- 
totics are governed by the (chiral) conformal dimension, 
A, of the operator appearing in the form factoid. In 
particular, if we write 

lim f^(6i,e 2 ) ~e y *\ 0i \, (4.30) 

\9i\^oc 

the quantity, y^, is bounded from above by A. This bound 
will prove useful in checking the end result. 

C. Basic Properties: One Particle Form Factors 

One particle form factors are in a sense trivial; Lorentz 
covariance completely determines their form. If ipfa) has 
Lorentz spin, s, then 

(0) = <O|V(O)4(0)|O) =ce s \ (4.31) 

where c is some constant. To determine c we use the 
theory's bound state structure. 

If a particle A3 is a bound state of two particles Ai 
and A 2 we can represent it formally as product of two 
Faddeev-Zamolodchikov operators: 

igl 2 A\(8) = res 5=0 A\(8 + 5 + iu%) A\(8 - iu\ % ). (4.32) 

Here res,5 = o denotes the residue at 5 = and T is the 
charge conjugate particle of 1. Implicit to this relation 
is the particle normalization (8\8') = 2n5(8 — 8'). The 
m's are given by u 2 ^ = ir — u 2 ^ where uf 2 is the location 
of the bound state pole of the S-matrix for particles 1, 2 
indicative of particle 3. g\ 2 is the amplitude to create 
particle 1 from particles 2 and 3. g\ 2 is defined by 

S \l{e) ~ i/^J-. (4.33) 
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Then we have 

^12/3" (#) = res <5=o ff 2 { 



iu\^9 + 5 + iu\y). (4.34) 



Thus knowledge of the two particle form factors com- 
pletely determines their one-particle counterparts. 



D. Two particle form factors for the SO(6) currents 

There are two possible two particle form factors for the 
50(6) currents, one involving two fermions, 



vfZZM) = (£#(0)4(02)4(0!)), 

and one involving two kinks, 



(4.35) 



$(81,82) = (G^ b (O)4(0 2 )4(0 1 )), (4.36) 

where the two kinks have opposite chirality. The con- 
tribution of the former to the spectral function is only 
felt at energies lu > 2 3 / 2 A while the contribution of the 
second is seen at u> > 2A. We thus focus on the latter. 

We begin by identifying the group theoretical structure 
of J a \: 



C^u82) = {Ca ab ) a$ U{8 1 ,6 2 ) 



(4.37) 



where C is the charge conjugation matrix introduced 
previously. (G<r ab ) a p is antisymmetric in a, b. Now 
{Ca ah ) a f } is not the only obvious choice of tensor. One 
also has (a ab C) a p, antisymmetric combinations built up 
out of j a Cj b or Cj a C-f b C, or some combination of all 
three (but not a ah as this choice violates obvious U(l) 
conservation). Perhaps the most natural choice is to 
make ^f^s explicitly symmetric in C: 

„/$(*!, fc) = {{Co ah ) aB + (<7 ab C) a p)f f t(0 1 ,6 2 ). (4.38) 

However this forces ju/^ -1 to zero, something we do 
not necessarily expect. We instead arrive at the choice as 
given above by explicitly checking the other possibilities 
for invariance under SO(6). Of the choices only (Ca ab ) a p 
is consistent as such. 

Having specified the coupling of a kink+anti-kink to 
G ab , we use C-symmetry to specify the corresponding 
anti-kink+kink form- factor ^,fgg(8i,6 2 )- Under the ac- 
tion of the charge conjugation matrix, C, the currents 
G" fc transform as 

CG ab C^ = {-l) a+b G ab . (4.39) 

We thus expect 

nf&VM = -(-ir+ b (-ir+& ^{e u e 2 ). (4.40) 

An ansatz for fj,f£g(8i, 2 ) satisfying this condition is 

„/$(*!, 82) = (O^W/^i, 8 2 ). (4.41) 



In particular we note the absence of any phase in the 
expression for ^f^(8i,8 2 ) relative to that of ^/^(fli, 9 2 ). 
As before Lorcntz invariance demands 

U(8i,8 2 ) = (e^+W 2 - (-lTe-^+^/ 2 )f(6 12 ). 

(4.42) 

Then by the scattering axiom, f{6\ 2 ) must satisfy 



3/2 - 6/iTT" 



G c (x) 



exp 

e x/ 2 



—G c (x) sinh(— ) 



f(8); 



sinh(x) 

The periodicity axiom takes the form 



(4.43) 



M/fa(02, 81) = -,J^i - 2™, 6 2 ). (4.44) 

The minus sign arises from Q a 0R a ,G ab = — 1- As shown 
in Appendix B, Q a p = — 1 and R a .c ab = 1 for neutral 
currents and vice-versa for currents carrying charge. 

In terms of the scalarized form factor, /, the periodic- 
ity axiom then reads 



f(-6) = f(6-2m). 



(4.45) 



Two kinks of opposite chirality do not form a bound 
state. Thus the form factor should have no poles in the 
physical strip. Thus with Eans. (|4.42|l . (|4.43|) . and l|4.45|l . 
the form factor is 

vfa b p(8i,9 2 ) = 

= A c {Ca ab ) a - p (e^ +e ^' 2 - (-l)M e -(*+«»)/a) 



x exp 



x s(x) 2tt 



(4.46) 



The hermiticity condition for this current form-factor 
reads 

f$j(8 u e 2 )* = ff a (6 2 - m, 6, - i*)C^C &a (4.47) 

Applying this constraint shows that A c is some purely 
imaginary constant. 



E. Two particle form factors for SO(6) symmetric 
tensor fields 



Here we will compute the two particle form factor cou- 
pling to an operator that transforms as a scalar under 
Lorentz transformations and as the 10 dimensional sym- 
metric representation under 15*0(6) (we indicate this op- 
erator as ip^ abc ^ A where a, b, c are anti-symmetrized). We 
also assume that the operator's conformal dimension is 
A = 3/8. Such an operator arises in the computation 
of the spin response for the ladders when considering 
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fermion bilinears of the form c Rb ^CLbi or antisymmetric 
combinations c Rb ^ c^ ab ^ — c Rab + Cn,i- Each of the fermions 
transforms under the 4 of 50(6) (5t/(4)). These bilin- 
ears are then in the 10 appearing in 4r <E> 4l. Here the 
elements of 4 R are (cr&j , c Rb i , c R .ab] , c Ra bi ) while the ele- 
ments of 4 L are (c[ H , c[ bp c[ abl , c\ ab] ). We can see, for 

example, that Cr^c^i lies in the 10 as in its bosonized 
form, 



J 

-Lb] 



4 +3 +c/>_ s ) _ gi| (01 -02+^3-04) 



it carries all three of the quantum numbers in 50(6). In 
considering tp^ abc ^ A , we rewrite the 50(6) fields appearing 
in the bosonization as 



D i±4>/2 



X2±«Xl, 



and so work in a Majorana like basis. In this basis we 
need to know how to construct the 10 out of the 6. From 
Ref. l46l we know the antisymmetric product of three 6's 
yields the 10. We thus denote this field as ifM bc } A where 
the A indicates anti-symmetrization. 

At the two particle level, two kinks of the same chirality 
will couple to such an operator. We will consider for the 
moment \jM bc ^ A coupling to two even chirality kinks. This 
then leads to form factors of the type 

^{OiM = (V^ 4(02)4,(01)). (4.48) 



50(6) covariance suggests that f a o 
form 



abc] A . 



h,9 2 ) take the 



Cl(0 7 [V 7 c ]) a /3 

+c 2 (7 to 7 V 1 )a/3 + c 3 (7 [ V0 7 c] ) Q/3 



+c 4 (7 [ V7 cl O) Q/3 )/(0i 2 ) 



(4.49) 



Here the f abc ^ A is a function of 9\ — 9 2 because ip^ abc ^ A 
is a Lorcntz scalar. 

We point out that the 50(6) tensors in the above ex- 
pression are symmetric in a and /?. To fix the values 
of the various coefficients we have two arguments avail- 
able to us. By checking explicitly for 50(6) covariance 
we see that all of the coefficients but c\ are zero. We 
can also fix the values of Ci relying on U(l) charge con- 
servation alone. Suppose \a) = (1/2, ±1/2, ±1/2) and 
|/3) = (1/2, ±1/2, ±1/2). Then by rewriting the Majo- 
rana fields as Dirac fields, one readily finds that 



f [16c] A _ _ • r [2bc] A 
J a(3 ~ L JaP ' 



(4.50) 



provided b, c = 3, 4, 5, or 6. As can be checked, this yields 
c 2 = C3 ~ C4 = 0. We thus obtain in the end 

fS c]A (0i,e 2 ) = {c^ 1 b 1 ^Uf(e 1 ^e 2 ). (4.51) 

where we have absorbed a constant into f(0). 



To determine f(9) we apply the scattering and peri- 
odicity axioms. The scattering axiom constrains f(9) as 
follows: 



/(-*) = (J«o(*) - \u 2 {6))f{6) 

T(3/4 + 6»/2i7r)r(l - 6 /2m) 
~ T(3/4 - 6/2m)T{l + 9/2in) 



f(9) 



exp 



f°° dx 

I — G s {x) sinh(.T0/i7r) 
Jo x 



G s (x) 



. 2e -3^/4 smh(x/4) 
sinh(x) 



The periodicity axioms reads 

f{-O) = Q a0 R a ^f{e-2m). 



f(9); 

(4.52) 
(4.53) 



Both phases, Q a and i? Q ^,, are equal to 1 (Appendix B). 
With this we arrive at the following expression for f{9) 



f(9) = sinh(-)exp 



dx G s (x) . 9 . x , ... 



(4.54) 

Although two kinks can form a bound state, they do so in 
the antisymmetric (vector) channel and so do not couple 
to ^ abc ^ A . Thus there is no need to add additional poles 
to the form factor. We then can put everything together 
leading to 



C c]A {9u9 2 ) = M^^CUshM-f) 



2/3 

x exp 



^sin 2 ( — (in + 6 12 )) 

x s(x) 2w 



(4.55) 



where A s is a real constant (as determined by hermitic- 

ity). 

The UV asymptotics of this form factor are easily 
found to be 



lim f 



abc]/ 



11111 J aB 

1 — *oo 



2( 



(4.56) 



Given that the conformal dimension of f^ abc ^ A is A = 3/8, 
we see that the UV behavior of this form factor saturates 
its upper bound. 

We now return to the situation of two odd kinks cou- 
pling to fl abc \ A . We thus want to consider the form fac- 
tors 



fg c]A (6 1 ,e 2 ) = (^ ab ^Ap(e 2 )A &{ei) 



(4.57) 



Under charge conjugation, the fields, f^ abc \ A , transform 
as 



j[abc] A 

We thus expect 

f l ^ c]A (9u9 2 ) = (- 



a-\-b-\-c j[abc]A 



-(-i) 



^a+b+c^_^a 2 +p2 j?labc] A 



(4.58) 



h,9 2 ) 
(4.59) 



19 



We may satisfy this constraint by choosing fj?s^ A {Qi, 82) 
to have the same form as /jS (01, 82), i.e. 

fS c]A (e 1 ,e 2 ) = (c 7 [a 7V ] ) a/ 5/(M- (4.60) 

Here f(0 12) is the same as before as all the other axioms 
governing the form factor remain unchanged. 



F. Two particle form factors for SO(6) vector fields 

In this section we are interested in computing the two 
particle form factors for an operator, tp a , transforming 
as both a Lorentz scalar and a vector under SO (6). For 
specificity we consider such an operator arising from the 
following symmetric combinations of fermion bilincars 
c i?i>T CLab i + c Rab'[ CLb l' This operator has conformal di- 
mension, A = 3/8. At the two particle level, two kinks of 
the same chirality will again couple to such an operator. 
We first consider the case of two kinks of even chirality 
coupling to this operator. 

We thus consider form factors, f£ a g(6i, 8 2 ), defined by 



fZa/iVM = (^(0)^(02)^(0!)). (4.61) 

Covariance under SO (8) suggests that f£ g takes the form 

ftoc^lM = [ci(7 a C) a 3+C2(C 7 a U}fe(ei,e 2 ), (4.62) 

where c\ and c 2 are constants. (That c\ and c 2 are not 
more generally independent functions of Q\ and 62 is eas- 
ily seen; the constraints Ean. (|4.22J) and Ean. (|4.23|l do 
not allow it.) c\ and c 2 can be fixed easily. Suppose 
I a) = (1/2, 1/2, 1/2) and \0) = (1/2, -1/2, -1/2). Then 
by rewriting the Majorana fields as Dirac fields, one read- 
ily finds that 



f 1 
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(4.63) 



This in turn forces c\ = 0. We then set c 2 = 1 as we are 
uninterested at this point in an overall normalization. 
As the operator is a Lorentz scalar, we must have 



1) = / e (0ia). 



(4.64) 



f e (8) is then constrained by the scattering axiom. Using 
the form of the kink-kink S-matrix for kinks of the same 
chirality together with the anti-symmetry of {C"f a ) a p in 
a and (3 we obtain 



/eH 



T(l - 8/2iir)T(8/2iir + 3/4) 
T(l + (9/2i7r)r(3/4 - 8/2in) 

1/4 + 6>/2z7T 



-1/4 + 0/2MT 



/«(«) 



exp 



— G v (x) sinh(— ) 



fe(0); 



<W = T3^(W- eX/2 ) 



-8x/2(l _ e 2x ) V (4.65) 



We now apply the periodicity axiom. It takes the form 
/eVM) = Q a pRa4,fLp(e - 2m). (4.66) 

As discussed in Appendix B, we find that for even kinks 
the product of phases QapRaip is equal to 1. With this 
the periodicity axiom reduces to 



U-0) = -/ e (0 - 2m). 



(4.67) 



The presence of a minus sign (compare Ean. (|4.58(l ) re- 
sults from the antisymmetry of C"/ a . The periodicity 
and scattering axioms then imply that f e (8) must have 
the minimal form 



f e (6) = sinh(0)exp 



dx G v (x) . 2/ x /■ 1 n\\ 
x sinnfi) 2tt 



(4.68) 

As two kinks of the same chirality form a fermionic bound 
state, /" o should have a pole at 82 = 0i 



in/2. Thus f" aB becomes, 

cosh(0i 2 J 



l a/3 



exp 



dx GJ 

r^VrSin 2 — "r + 012 

x smn(i) Zir 



,(4.69) 



Here A v is the some normalization with mass dimension 
[m] 3 / 8 . We now consider the situation of two odd kinks 
coupling to ip a . We thus want to consider form factors 
of the type 



/^(*i.*a) = «A0(0 3 )^ (fll) ). 



(4.70) 



All but the periodicity axiom is of same form as for the 
two even kinks. For the two odd kinks, the periodicity 
axiom changes to 



f o (-0) = fo(0-2m). 



(4.71) 



The difference with Ean. H4.67fl arises because the OPE of 
the odd kinks with the operator produces an additional 
sign. This sign changes the form factor involving two odd 
kinks to be 



-4, 



exp 



sin(7r/4) 
dx G v (x) 



sinh(0 12 /2) 
(C7 ) a g . ,. . X 



cosh(0i 2 ) 



sin 2 ( — (mt + 12 )) 
x sinn(i) Z7r 



, (4.72) 



Here A v is the same normalization as for the even kinks. 
We can fix the normalization of /° a s(^i>^2) relative to 
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fea/3 (^i ' ^2) because both form factors must have an iden- 
tical bound state pole structure. This then mandates the 
relative factor, To determine the phase of A v we 

again employ the hermiticity condition: 

fLpQuW = C m C &a f a e - pa {6 2 - in, 8, - in). (4.73) 

This implies A v is real. 

We can compute the asymptotics of these form factors. 
We find that for the even kinks, limg^oo /"^(^l, #2) be- 
haves as exp(36>i/8). Thus the bound is saturated as we 
found in the case of two kinks coupling to a field trans- 
forming as the SO(6) 10 dimensional symmetric repre- 
sentation. However for the odd kinks we find instead 
that lims^oo f a oaB {8i,e 2 ) = exp(-l/80i). 

If we had considered an operator different than 



c -Rb'\ CLab i + c RafcT Cifc i' still transforming as a SO(6) 
vector, we would arrive at the same results with the pos- 
sible caveat that the form factors of the even and odd 
kinks may be swapped. For practical purposes this is a 
distinction without a difference. In any sum over form 
factors that would appear in a computation of a correla- 
tion function, the form factors involving both even and 
odd kinks would appear with equal weight. 



G. One particle form factor of 50(6) vector fields 

The SO(6) vector field will couple to a single Gross- 
Neveu fermion, a. So we consider the form factor, 

tf(9) = (^(0)4(6)), (4.74) 

which must take the form 

h a (9) = cS a b . (4.75) 

To determine the constant, we use the two particle 
form factor f« p . Let \a) = (1/2,1/2,1/2) and \/3) = 
(1/2, —1/2, —1/2). A a , a = 1, can be written in terms of 
these two even chirality kinks: 

ig^A a (e)=res s=o A a (0 + 5 + iu% a )A a (e-iui a ). (4.76) 

Again the u's mark out poles in the S-matrix indicative 
of bound states. Here they are given by 

u$ a = €a = K- i4 = tt/4. (4.77) 

g^p can be determined up to a phase from the kink-kink 
S-matrix as before to be 

/2 r(7/4)\V 2 

•^U^wiyJ ■ (478) 

Given that, 

iS&fifW) = «*,=«, ±f^(e-in/4,6 + S + in/4), (4.79) 



together with the hermiticity constraint, 

= (4-80) 
we find (up to a sign) the constant in Eon. 1(4. 75(1 to be 



y a f3 



exp 



dx G*^(x) . ^2/-^ 
x sinh(a;) 4 ' 



.(4.81) 



If we repeat the above procedure using the two particle 
form-factor involving t 
obtain the same result. 



form- factor involving two odd chirality kinks, we 



H. Summary of form factor results 

In this section we summarize the results of this section 
for the various form factors. 

Two particle form factors: 

The form factors of the SO(6) currents, G° b , are: 



= A c {Ca ab ) a - p (e^ +0 ^' 2 - (-1)^-^+02)72) 



x exp 



dx G c (3/ ) .0/3//. ^ \\ 

rv sm (5- i7r + 0i2 

x s(x) In 



(4.82) 



The form factors of SO (6) symmetric tensor fields, 

^[abc]A ; are . 



x exp 



dx G c \X ) o / 33 / „ xx 

iLisin 2 — (m + 9 12 )) 

x s[x) In 



2 

(4.83) 



The form factors of SO (6) vector fields, ip a , involving 
two even kinks are 



exp 



cosh(6»i 2 ) 

dx GJ 3/1 .9/33. . . 

T^ Bm ^ tv + ^ 



,(4.84) 



and for two odd kinks change to 



sinh(0 12 /2) 



p sin(7r/4) cosh(6'i2) 



exp 



dx GJ 

riV4sin 2 (— (m + 12 )) 

x sinn(i) In 



(4.85) 
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We note that depending on the particular vector field 
considered the form factors of even and odd kinks may 
be swapped. 

One particle form factors: 

The one particle form factor for SO (6) vector fields is 



where c is given by 
• \ 



9*p 



■ exp 



Sbifi) = cSl (4.86) 



-^P- S inli 2 (-) .(4.87) 
o x smh(s) 4 
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APPENDIX A: CALCULATING THE 
DYNAMICAL SPIN-SPIN CORRELATORS 

In this appendix we present the detailed derivations of 
the functions Ji,2,3(w, k) in Eans. (|3.1U|) - (|3.12|) that de- 
termine the dynamical susceptibility In subsection A 4 
we estimate the amplitudes Aij that determine the rela- 
tive spectral weights of the various contributions to the 
structure factor. 



V. SUMMARY AND CONCLUSIONS 

In this work we have calculated the dynamical mag- 
netic susceptibilities of doped Hubbard-likc ladders in 
the framework of a low-energy description by means of 
the SO (6) Gross- Neveu model. The most prominent fea- 
ture of the low-energy spin response is a narrow peak at 
wave-vectors (along the leg direction) Kp a b + Kpi> and 
—Kpab — Kpb. This peak corresponds to a fermionic 
bound state of two SO (6) kinks, which is broadened by 
the gapless bonding charge mode. This narrow peak 
sits on top of an incoherent scattering continuum of two 
SO(6) kinks. The low-energy dynamic magnetic response 
in the vicinity of the wave numbers ±2Kpi, and ±2Kp a b 
is rather featureless. The response at low momenta and 
in the vicinity of ±(Kp a i,— Kpi,) exhibits incoherent scat- 
tering continua with a threshold singularity. 

Other dynamical correlation functions such as the 
single-particle Green's function, density-density or super- 
conducting correlators can be determined along the same 
lines. 
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1. Calculation of Ji(w,fc) 

The contribution Ji(w,fc) involves fields arising from 
right /left symmetric combinations of fermion bilincars of 

the form c^^cr^i + c ihj c L&i- These fields in turn are 
expressible purely in terms of zeroth Lorentz component 
of the the SO (6) Gross-Neveu currents. For example 



V 2 r^3R + ^Il^sl 



-(—iGl 6 + G 



J 45 



G 



36 



*G° 5 ). (Al) 



We can then write down J± as 



Jx(u;,k)= / drdxe*-™ (TG^(x, r)G° d (0)>. 



(A2) 

Here cd gives one SO(6) current. Given that the differ- 
ent currents that arise in expanding c^cr^ + c^^CLbi 
give identical contributions (up to a constant) we simply 
write J\ in terms of a single current, understanding that 
any degeneracy is fixed through absorption into the as- 
sociated pre-factor, Aji, j = 1, 2,3. Under a form factor 
expansion, the lowest energy excitations contributing to 
Ji are two kinks of opposite chirality. We first focus on 
the imaginary piece of J\ . For the real piece it is easiest 
to obtain via a Kramers-Kronig transformation. 
The imaginary piece takes the form, 



/dO dO 
±±*\ffi(6 x ,e i )\ i 8{u,-E<fi u o 2 )) d(k-p(e u e 2 )) 



(A3) 



where R/Lf%fa(6i> $2) is the kink-kink form factor for the current operator computed in Section 4 and 

E(0 1 ,Q 2 ) = m [cosh 6»i + cosh 2 ] , 
P(8 U 6 2 ) = — [sinh X + sinh 9 2 ] . 

Vp 



(A4) 
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FIG. 5: The real and imaginary parts of the function 
Ji(Lo,m/vF) for low frequencies. For general values of q there 
is a square root singularity in the imaginary part above the 
threshold, ui = y v 2 F q 2 + 4m 2 . 



We mean the indices, a and /?, of the sum, Yla0i ^° run 
over both kinks and anti-kinks. Performing the sum, we 
obtain 



ImJi(LU,k) = 87r 



I d9 x d9 2 , f0/Q 
(2tt) 2 



x<5( w - E(9 1 ,9 2 )) S(k - P(6 1 ,6 2 )),(A5) 

where (up to a constant which will be absorbed into the 
appropriate An), 



\f°(0i,0 2 )\ 2 = Sinn 
dx G c (x) 



2/9\ +02-. 

' 2 ' 



x exp 



(l-c(a;)cos( ) 

p X S(X) IT 



(A6) 



and G c (x) is given in Ean. (|4.46|) . We then perform the 
integrations over 9\ and 02 obtaining 



ImJi(w, k) 



8v F k 2 6{u - Vfc 2 +4m 2 ) 



( w a _ fc 2)3/2 ( w 2 _ fe 2 _ 4m 2)l/2 



x exp 



dxG c (x) ,0i 2 x 

rfi-w cos ( — ) 

X S(X) IT 



where 9\2 is given by 



?i2 = cosh 1 (- 



2m 2 



2m 2 



■)> 



(A7) 



(A8) 



and k = vpk. We plot ImJi(o;, m/vp) in Fig0 



2. Calculation of J2(u),k) 

We now consider the computation of J 2 . This contri- 
bution arises from considering fcrmion bilincars of the 
form c^c/^. These bilinears are expressible as prod- 
ucts of the gapless charge degrees of freedom with fields 



transforming under a 10 dimensional symmetric repre- 
sentation of SO(6). For example 



be]. 



(A9) 



ibe 



Here the fields i{j^ abc ^ A are denned and discussed both in 
Sections 2 and 4. The coefficients c a i, c do not need to be 
determined exactly as the contribution of each ifM bc ^ A to 
J2(w, k) will be identical up to some constant. Any unde- 
termined constant will be fixed in the end by appealing 
to a sum rule. 

J2(x,t) is then the product of two correlation func- 
tions, one coming from the gapless charge correlations 
and one coming from the gapped SO(6) sector. We write 



where 



J2{x,t) = J 2g {x,T)J c (x,T), 

J c {x,t) = (T e -30c + (z,r) e ^ c+ (o,o)^ 
= A c {x 2 +t 2 ) k ' a - 
J2 9 (x,t) = (T4,^ A (x : r)^ abc ^(0)). 



(A10) 



(All) 



Here K is the Luttinger parameter of the gapless charge 
excitations and A c is some constant which we will absorb 
where appropriate into the ^4's of Eans. (|3.1UI3.12|) . We 
evaluate J2 g in the same fashion as J\, namely we use 
employ a form factor expansion that is truncated at the 
two kink level. We thus obtain 



J2 9 (x,t) 



In 



a/3 



d9 1 d9 2 , ,abc A ,n a \\2 
T a p (#1, 02 J I 

z ,iT{-iu+t)-\T\E{0 1 ,8 2 ) p ix{k- S i g n(T)P(8 1 ,e 2 )) 



A£/d0-l/^ CA (20-)l : 

aj3 



x A" o (2mcosh(0_)((2;/u F ) 2 + t 2 ) 1/2 ). (A12) 

Here a and (3 are kinks of the same chirality. In the 
second line of the equation we have a made a change 
of variables to 6± = (0i ± 2 )/2 and performed the 9 + 
integral. Using the form factors of Section 4.4, we can 
write |/^(0i, 2 )| 2 as 



x exp 



2)| 2 = sinh(f f 

^%W (1 _ c(x)cos( ^: )) 

q X S[X) ' 7T 



(A13) 



Substituting Ean. ljA"T3|) and Eqn.JXTU into Ean .lETTl . 
we can carry out the x — t integrals following Refsl52land 
IFlI Up to some constant (which we absorb into A12/A22), 
we obtain 
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J 2 (w,fc) = 



v F 



K 



sinh 2 (6>) , 

ae 1- 

■ (cosh6<)~ 4 



A' 

~ T' 1 ' 



- esign(w) — k 



exp 



dxGJx), . . ,20x„ 
— 1-cx cos 

X S{X) IT 



(A14) 



where 




2 3 4 
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3 (x,t) = c f Rab1 -(x,T)c Lbl (x,T) + ab . (A16) 



The above symmetric combination of fermions is equiva- 
lent to the antisymmetric field tp a discussed in Section 4. 
This becomes clear under bosonization of these fermions 



<D 3 (x, t) = 2e~i^ + ^ cos(i(0 3 + 4 )). (A17) 



FIG. 6: The real and imaginary part of the function, Ji(u), 0), 
for low frequencies and K = 0.95. 

We plot the imaginary part of J 2 (w,0) in Fig|S|for the 
special case K = 1. 

3. Calculation of Js(u},k) 

Finally we consider the calculation of J 3 (ui,k). J 3 
comes about from considering the correlations of fermion 
bilinears of the form Cr^cl^ + Cjj b ^Cx, ^. In particular 
we have 

J 3 (w,fc) = J dTdxe^+^ kx (O 3 (x,T)Ol(0,0)) , 

(A15) 



The model is symmetric under the Z 2 transformation 
6*3,4 — * —6*3.4. It is clear that this field couples to the the 
Gross- Neveu fermion, exp(i</) fi / i2 ), whereas the antisym- 
metric combination of Fermi fields, c Rab ^CLbl — c mi CL aH' 
proportional to sin(i(6 ) 3 + #4)), does not. The 6 is then 
identified with symmetric combinations of fermion bilin- 
ears. 

J3 has a structure similar to that of J 2 : it is a product 
of correlation functions of both the gapped and gapless 
degrees of freedom, 

M*,t) = J 3c j{x,t')J c {x,t). (A18) 

However the gapped part of the product is different in 
that it has contributions at low energies from both single 
and two particle form-factors: 



J 3g (x,r) = ]T|c a | 2 / |i|^ O(0) |4 (ei)) |^-kl^c( 9l ) + ^ ig nM^ s(ei ) 

a J n 

^E^E / ^l«(0)l^(^K(^))| 2 e-M^ 



e 2 )+ixsign(T)P(e 1 ,d 2 ) 

a/3 

J 3 g,l(x,T) + J 3g , 2 (x,T) . (A19) 



The contribution of the one-particle form-factor marks spectral function. We do note that this modes broadens 
the most significant difference between J 3 and J 2 . It in- out because of the gapless charge excitations. The contri- 
dicatcs a propagating coherent mode contributes to the bution of this mode has been previously studied in Ref. 
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l25l We now drop the sum. ^ a |c a | 2 , as each term gives 
an identical contribution up to some constant. Using the 
form factors of Section 4.5 we can then write 

hg,i{x,T) = A 2 A -K (V2m{T 2 + {x/v F ) 2 )) 



1 f 

d = — exp[-2 / 
9 Jo 



dx G v (x) 
o x s(x) 



,tanh(26L) 2 sinh(0_ 



J 3s , 2 (*,r) = A\- 2 jd9_{ , cosh(2 0_ )2 . 

xK (2mcosh(9-)((x/v F ) 2 +t 2 ) 1/2 ) 



x exp 



dx G v (x) 29-X 
ryl-cW cos 

X S{X) 7T 



(A20) 



Here the relative weights of J 3ffj i and J 39! 2 are fixed 
using the bootstrap as described in Section 4.6. g is given 
by 

2= VHX7/4) 
9 3r(5/4) ' 

via Eon. (|4. 78(1 . Again following Refslfi2l and lfill we can 

rewrite J 3 (o;, fc), up to some overall constant, in the form 



J 3 (w, k) 



v F 



2g 2 m 2 
2 K '*v F 



exp 



dx G v (x) ,o,x.\ ( K K uj 2 — k 

^Vrsinh 2 - mi ,1 ,1, ^ 

x smh{x) 47 V 4 4 ' ' 2m 2 



dO- 



1 



(cosh (9)- 



tanh(26») 2 sinh(6») 5 



cosh(2<9) 2 
x exp 



/ _ K _K to 2 -k 2 
1 ^ 4' 4' ' 4m 2 c 2 (0) 
3 dxGJx) , , , ,2Qx s ' 
a; s[x) 7r 



4. Estimating the Amplitudes Ai 



(A21) 




2 3 4 
co/m 



FIG. 7: The real and imaginary part of the function J3(w,0) 
for low frequencies and K = 0.95. There is a singularity 
above the threshold at lj = \/2m. The kink-kink continuum 
contributes for u > 2m. 



We plot J 3 (u), 0) in Fig[7|for the case K = 1. The domi- 
nant feature is a narrow peak corresponding to the SO (6) 
fermionic bound state. The peak is not sharp due to the 
admixture of excitations from the gapless charge sector. 



We now turn to estimating the various normalization 
constants, i.e. the A's appearing in Eqns I3.1UI3.121 In 
the present context, perhaps the best way of determin- 
ing them would be to compare the held theory results to 
numerical computations following the lines of BSi. How- 
ever, such computations are presently not available. A 
crude way of determining the order of magnitude of the 
A's is to compare our results to the susceptibility calcu- 
lated in a Random Phase Approximation. The latter is 
given by 



X 



RPA 



(w,q) 



Xab{uJ,q x ) 



+2cos(q'j / )xi2(w,fe) 

.-(0) 



l + 2U X ^(^,q x 



(A22) 



a b 



Here x ( w i Qx) is the susceptibility matrix of non- 
interacting electrons (U = 0) on a two-leg ladder. In 
order to fix the constants Ay, we now require 



du , 



Im X ^(c,q) 



dio 



Imx(w,q), (A23) 



where A is an appropriately chosen energy scale. The 
rationale behind Ean. l|A23|) is the following. We don't 
expect the RPA to give an accurate description of the 
susceptibilities at low energies as it for example fails to 
account for the dynamical generation of spectral gaps. 
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However, we expect RPA to give a reasonable account 
of the static susceptibility (which corresponds to taking 
A — > oo ) and hence also of the integrated spectral weight 
at low frequencies. Field theory does not give a good 
account of the static susceptibility, which is reflected in 
the fact that one cannot take A — ► oo in the field theory 
calculation as the resulting integral would diverge. We 
therefore arbitrarily set A = 10m. A useful consistency 
check is that the results for Aj are essentially unaffected 
by small changes in A. We choose to consider Ean. l|A23jl 
rather than the integrated spectral weights in order to 
suppress as much as possible the neglected contributions 
that occur at higher energies in the field theory calcula- 
tion. On the field theory side we have 



doping 


m SO{6) (Wi/2filling) 


K 


1% 


.47 


.94 


5% 


.29 


.915 


10% 


.21 


.93 


25% 


.13 


.945 



In this table we have additionally listed the correspond- 
ing Luttinger parameter for the different dopings, also 
available from Ref. l|43T) . Having fixed vp and m in the 
low-energy effective SO(6) theory in terms of t and U 
we now can use the comparison between the field theory 
and RPA calculations of the integrated susceptibilities to 
estimate the amplitudes Ay. We find 



1 T , n 3{An+A 21 ) 
— Imx(w, q, -) ps Ji(w, q), 

TT Z 47T 

-ilm X (c,2A Ffc ,J)«^J 2 (c,0), 
tt 2 An 

-~hnx(u, 2K Fab , ^) ps J 2 (cj, 0), 

TT Z 47T 

-Imx(w, K Fb - K Fah + q,^r)^ ^^-Ji((J, q), 
TT 2 An 

--lm X (u,K Fb + K Fab ,?-)tt ^ij 3 ( w ,0). (A24) 

TT 2 47T 



The relevant frequency integrals arc 



doping 


An 4 


- An 


Al2 


A22 


A31 


A32 


1% 


0.24 




0.13 


0.10 


0.14 


6.61 


5% 


0.25 




0.045 


0.070 


0.14 


4.70 


10% 


0.26 




0.12 


0.091 


0.22 


3.24 


25% 


0.61 




0.091 


0.79 


0.24 


4.94 



A table of the field theoretic amplitudes, Ay, 
computed for U = t and t± = 1.5t. 

All values for the As are in units of 

We thus have all the ingredients at hand to calculate 
the spin-spin correlation functions of the doped ladders. 



f 10m dw 5m ^ v F 
/ — Ji (w, ps 1.04—, 

f Wm dcu . 5m x v F 
/ — J 2 (w,— ) ps 1.79-1, 

f 10m dw T . 5m N n „ « f , 4 „ rN 

/ — J 3 (w, — ) ps 0.23^. (A25) 

In order to compare to compare to the RPA calculation 
we need to estimate the gap m as a function of t, t± and 
U. We do this in two steps. Firstly, the spin gap at half- 
filling A s (l) has bee n deter mined numerically for small 
values of U/t in Refs |2ll24j . For t± = 1M and U = t it 
was found that 

A s (l) ss 0.05i . (A26) 

In order to determine the kink gap m we need to know 
how A s evolves under doping. As numerical results are 
not available in the literature for small U, we instead 
use the results of Ref. lf43T) . where the doping depen- 
dence of the spin gap was analyzed within the frame- 
work of a low-energy description of the undoped ladder 
in terms of the SO (8) Gross-Neveu model. At half- filling 
we have mso(8) = A a (l). Assuming that the Fermi ve- 
locity is well approximated by the non-interacting value 
v F ps 1.32t (for t± = 1.5t), we find that the gap m 
for SO (6) kinks at four different doping levels is roughly 
equal to 



APPENDIX B: DERIVATION OF FORM FACTOR 
AXIOMS FOR GENERAL OPERATORS IN SO(6) 
GROSS-NEVEU 

In this appendix we derive the form-factor axioms for 
SO (6) Gross-Neveu for general semi-local operators. We 
take the following approach. We first consider the gen- 
eral form-factor axioms as written by Smirnov in Ref. 
(|40|) for the current and stress energy operators. These 
are written for a particle basis which is not C-symmetric. 
In this basis the form-factor axioms are not easily gen- 
eralized. To surmount this difficulty wc transform to a 
particle basis which is C-symmetric. Here the form-factor 
axioms can be readily written down for all semi-local op- 
erators. In particular the semi-locality phases appearing 
in the axioms can be computed using operator product 
expansions. We then recover the form- factor axioms for 
general semi-local operators in the original basis by in- 
verting the transformation. As an important consistency 
check of both the transformation and the operator prod- 
uct expansions, we recover the axioms for the currents 
and stress-energy tensor as originally stated. 

1. The N-Particle Form Factor Axioms for the 
Currents and Stress-energy Tensor in SO (6) Gross 
Neveu 

The N-particlc form factor axioms for the currents and 
the stress-energy tensor in SO (6) Gross-Neveu are given 



2G 



by Smirno\isi as follows: 



f°(. . .,9 j+1 ,6j, 



/"(fll, . . . , 0„ + 27rf) 6l ,.„, eB = (-l)^J l[ei] f°(6 n ,6 U . . . , fln-l)e„, 61 ,... ie „_i; 

iRes| en=en _ i+i7r / (^,...^ n ) ei ,... >en = f°(6 u ...,e n - 2 ) el> .. Wn _ 2 C en ^ t 

x '4 • • • S^zl - (fln _ a _ _ 02) ... 5 £ 7 ;- 3 el- 2 (0„_ a - 0„_ 2 ) . (Bl) 



The first equation is the generalization of Eon. (|4. 2211 of 
scattering axiom to matrix elements with N-particles. 
The second equation is the N-particle generalization of 
Ean. (|4.23fl . the periodicity axiom. The phase that ap- 
pears on the r.h.s. of the periodicity axiom is defined 
through 

!1 for even parity kinks 
2 for fermions (B2) 
3 for odd parity kinks, 

where the l[ej] are the rank of the representation the 
particle, ej, falls in. For the operators to which these 
axioms apply, the phase is always either ±1 as X^[ e j] 
is always a multiple of 4. While the above axioms are 
completely general, we will concentrate solely on form 
factors involving particles which are either kinks or anti- 
kinks. The form factors have simple poles at the points 
d n = 9j+iir. The residues of these poles are the subject of 
the final axiom. Here C is the charge conjugation matrix 
and by necessity l[e n ] +l[e n -i] = 4. The annihilation pole 
condition relates form factors with different numbers of 
particles. In the annihilation pole axiom a phase similar 
to that of the periodicity axiom appears. 

From one viewpoint, the phases in the latter two ax- 
ioms can be justified by insisting the currents and the 
stress-energy tensor satisfy both certain properties un- 
der charge conjugation and appropriate commutation re- 
lations. From another, the phases can be seen as a result 
of the fact that the C is not symmetric, i.e. C 2 = — I, 
and the non-trivial semi-locality between the non-Cartan 
currents and the fundamental kink fields. Understanding 
the phases in this fashion provide the most ready path 
to generalize the axioms to other fields. But because the 
phases have two sources, this path needs to be indirect. 
The first step is to transform the particle basis to a ba- 
sis where the corresponding charge conjugation matrix is 
symmetric. In this basis the form-factor axioms can be 
written down readily - the only additional phases that 
appear are semi-locality factors and these are well under- 
stood. To obtain the form-factor axioms in the original 
basis we then invert the transform. 



2. Transformation to a C-symmetric basis 

The path we then intend to follow is the one Ref. 
uses for SU(2). There the transformation to a C- 
symmctric basis is straightforward because of the simpli- 
fied particle content. Here, in contrast, the transforma- 
tion is more involved. We define the transformation to a 
C-symmetric basis as follows. A n-particle state behaves 
under the transformation, T, via 

i ei ---I e „|0) = TA ei ---A en \0) 

= (l[P eiej )A €l -..A en \0), (B3) 

i<j 

where the P^i = ±1 are given by 

Pa/3 = P/3& = Pap, a P\ 

p P0 = - p pp = -1 ; 
p — — p- - — i • 

Pafi = -P«fK (B4) 

and 

PapPpa = -1 if a ^ /?; 

P & pP(3 & = -i if « + 0. 
To meet this last set of conditions, we arbitrarily set 

P+++, — h- = P+++,+ = h = P-+-,-\ 

= P+— ,—+=P— +,-+-= -1. 

Here a, (3 correspond to even chirality kinks while a,/3 
correspond to odd chirality kinks. This choice of the 
transformation is not unique but is equivalent in effect 
to other possible choices. 

In moving to this transformed basis, the S-matrix also 
undergoes transformation. The transformed S must sat- 
isfy 

A ei A e2 = SeleZ A e > 2 A^ (B5) 
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This implies the relation 



we require that for any £3 that 



C e l e 2 Q € l € 2 p r> , , 

J £i£2 — J £l£2 - r eie 2 Jr e 2 e 1 



(B6) 



By construction S is C-symmetric. In particular it satis- 
fies 



- *) = S^S(9)C^C aS , (B7) 

where C 1 = —47274761 or equivalently, C a& = C aSl P a a, is 
now a symmetric matrix. 

This transformation, as constructed, satisfies an im- 
portant consistency condition. Considering the relation 
involving three particle states, 



A ei A e2 A e3 — SelezA^A^A^, 



(B8) 



we sec that to transform this relationship covariantly to 
the original basis, i.e. 



AAA - C e i e 2 AAA 



(B9) 



P P — Pi Pi 



(BIO) 



holds for any ei, £2, ei, e 2 such that S\\\ 2 is non-zero. The 
transformation as given above does indeed satisfy this 
condition. This condition guarantees both that S satis- 
fies the Yang-Baxter equation and as we will see, that we 
can recover the axioms given in the original basis from 
the axioms stated in the C-symmetric basis. 



3. Form Factor Axioms in C-Symmetric Basis 



In the C-symmetric basis the axioms for arbitrary 
semi-local operators have the following form: 



iRcs 



f c 

-1+27T J 



' * " ' j + 1 ' j 1 " ' 

+ 27ri) ei ,..., e „ 

■ ■ ■ j 9 n )ei ,...,€„ 



s 



7+1 ) / (• ■ • >0J)03+1) ■ ■ •)■ 



-r.; f u (0 nj 9±,..., ^n-l)e„ 1 ei,...,e„_i ! 



2iriw(0,e„) fOf 

/ (#1> • • • j 0n-2)ei,...,e^_ 2 C , £„,e^_ 1 

n-1 — $l) -§72 lii^n-l — ^2) ' ' ' Se^-itn-z (#n-l 



?n-2 y 



(Bll) 



The phase factor e 27 ™ w ( e, > < arises from the semi- 
locality of the fields in the SO(6) Gross-Neveu model. 
Working in Euclidean space, fields can be said to be semi- 
local if their product sees a phase change when the fields 
are taken around one another in the plane via the an- 
alytic continuation z —> ze t27T ,z — > ze _l27r (in conven- 
tions where right-moving fields depend only on z) . More 
specifically we have 

V> e (ze l27r , ze- l2 *)O(0) = e 2 ™^°' e ^ e (z, 1)0(0) (B12) 

Here tf) e is the field that is associated with the particle 
A\ (9) . ip e is chosen in such a fashion that it is both self- 
local and bosonic. For the case at hand we represent the 
kinks through the vertex operators 

Arun.n «-» ei'- 71 * 2 ^ 2 * 3 ^ 3 ^', (B13) 

where here a = ±1 and <fii, i = 2, 3, 4, are the bosons as 
discussed in Section 2. The choice of field, tp e , influences 
the phase, e 27riW ( c '^) j as different choices will have differ- 
ent operator product expansions (OPEs) with O. Our 
choice passes an important check: the phase, e 27riW ( c '^) j 



so obtained from this assignment allows one to recover 
the form-factors axioms as stated in the original basis. 

We now write down the semi-locality phases for the 
currents and the stress energy tensor. The stress en- 
ergy tensor, T^, is both chargeless and local. Its OPE 
with the kink fields as described above uniformly loads 
to e 27r M r>, V) = 1, The phases for the currents can be 
written in a compact form. Any current, J a( g can (group 
theoretically) be thought of as a product of an even and 
odd kink, a and {3. If the e corresponds to a kink 7 (of ar- 
bitrary chirality) the OPE then leads to the phase being 
given by 



e 2^(J„ ;5 , 7 ) = -p^p^p up (B14) 



This can be checked directly. With the phase in this 
form, we recover the form factors in their original form 
more readily. 
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4. Recovery of Form-Factor Axioms in 
Non-C-Symmetric Basis 

We now show that we can recover the form-factor ax- 
ioms for the currents and stress-energy tensor as given 
in the original basis (i.e. Ean. i|Bl|l ') by transforming 
the form-factor axioms in the C-symmetric basis (i.e. 
Ean. (|Bll|) '). This shows both that the transformation 
and our choice of fundamental fields to represent the 
kinks in the C-symmetric basis are consistent. 

The demonstration that we recover the scattering ax- 
iom, the first equation in Eqn.(rcfeBi), is trivial. It fol- 
lows immediately from the complementary fashion by 
which the S-matrix transforms relative to the basis (i.e. 
Ean. (|B6|) ) and the fact that the condition in Ean. (|B10H 
holds. 

The periodicity axiom provides the first non-trivial 
check of the transformation. We will first consider the 
topologically neutral operators, for example, the Cartan 
currents and the stress-energy tensor. In this case the 
semi- locality phase is 1. If we apply the transformation 
to the second line of Ean. (|Bll() . we see that to recover 
the periodicity axiom as stated in its original form, wc 
require that 



n 



(B15) 



To establish this we first note the any set of particles 
coupling to a topologically neutral operator can be de- 
composed (not uniquely) into the following three subsets: 

1. n e complete sets of the four even kinks. 

{(+ + +),(+- -),(- + -),(- -+)}; 

2. n a complete sets of the four odd kinks, 

{(- -),(- + +),(+-+),(+ + -)}; 

3. n p pairs of particle-antiparticles {a, <5}.(B16) 

We thus decompose the n-states as follows, n = An e + 
4n Q + 2n p . It is straightforward to show that the P's of 
the transformation satisfy the following conditions 



QEeven kinks 



QtEodd kinks 

PajP^aPa^P^a 1; v^^-'O 

where 7 is a kink of arbitrary chirality. The two sides of 



Ecin. (|B15(l can then be rewritten as 

?i-l 

Y[P enU P uen = +«P; 

These two are equal and we thus see that we recover the 
periodicity axiom for topologically neutral operators. 

The case of the non-topologically neutral operators, 
specifically the non-Cartan currents, is dealt with simi- 
larly. Denote by J a p the particular current of concern. 
Any n-particle state coupling to it will consist of some 
topologically neutral combination of n e complete sets 
of even kinks, n complete sets of odd kinks, and n p 
particle-antiparticle pairs plus one of the following: 

1. the kink pair a and /3; 

2. four even chirality kinks of the form 

{/3,/3 and two even kinks 7 7^ 7', 7,7' 7^ a,/3)}; 

3. four odd chirality kinks of the form 

{a, a and two odd kinks, 7 7^ 7 , 7, 7' 7^ a, /3}; 

(B19) 

Each of the three combinations listed above carry the 
quantum numbers of J a p . In case 1, n = 2 + 4n e + 4n D + 
2n p . To recover the periodicity axiom in its initial form 
we must satisfy 



n-1 



S 2*MJ«*.*») Y[ Pe n e t Pe z e n = (-l)*^' 



(B20) 



Using the expression for the semi-locality phase in 
Eqn. (|B14|I we see that 



2-7riw(J cl/ 3,e„) 



I Pene.Pe, 



while 



(_l)iE;jZkj] _ r^\n e +3n +n p +l ^ (B22) 



We thus see that in this case we recover Smirnov's ax- 
iom for the non-Cartan currents. The other two cases in 
(|B19j) proceed identically. 

The final axiom to consider is the annihilation pole 
axiom. If we transform the final equation of (|B11|) wc 
obtain 



)i-2 

iRes\ en=en _ i+i7r <>„) = / O (0i,...,^-2) ei ,...,<_ 2 ^11^^ 

i=l 
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n-2 



71 — 3 



X *S* 7l ei 



72^2 l^n-l — C2 



) -s?::?t 



>n-\ — Vn-2 



(B23) 



We first demonstrate that 

n-2 

Y[ P en6i P en _ iei = 1 (B24) 

i=l 

Given that e„ and e„_i are particle-antiparticle, the n — 2 
states remaining must be of the type given in Eqn. (|B16|) 
if the operator carries zero topological charge or of the 
form given in Ean. (|B19|) if the operator is a non-Cartan 
current. In either case this product over ti can be broken 
down into pairs of cither two even chirality kinks, {/3, 7}, 
two odd chirality kinks, {/?, 7}, or one even kink and one 
odd kink, {/3, 7} (where /3 is not necessarily the anti- 
particle of 7) . That the above product is 1 then reduces 
to the following equalities which can easily be checked: 

PapPaftPayPaj — lj 

PapPapParyP&i = (B25) 

where (3 and 7 are arbitrary even kinks. Noting that 
f-' n _ 1 = e n -i, we now demonstrate 

n-2 

i>j 
n-3 

X ( Pli^iP(-' i li-l)P^n-l ! in.-lPe.' n _ 2 1 n -3 

i=2 
n-2 

= II Pen-^Pe^ (B26) 
i=l 

If we can do so, the axiom reverts to the form in Ean. (|Bl|) 
if we can show 

n-2 

e 2^(0^ l} -Q p^ ^^ = (— l)iEr=i 2; [ e d. 
i=l 

(B27) 

But this equality is almost the same as Ean. l|B"T7|) and 
Ean. (|B22|) and can be established on identical lines. 

To show Eqn.{B26| holds, we appeal to the behavior of 
the transformation under scattering. The phase on the 
l.h.s. of Ean. i|B26() arises from the scattering (note that 
we consider only the case e' n _ 1 = e„_i) 

A en _ 1 ■ ■ ■ A ei = 5 7 "ri (Q n -\ — Oi)S l2 e% (0 n -i — 62) ■■■ 



x ~Sl::li: 2 2 (0„-i - 0n-2)A < 2 ■ ■ ■ A e ' i A e ^_ i . (B28) 
Applying the transformation, wc know this must become 

A en _ 1 ■ ■ ■ A ei = — Qi)S?(lll(6 n -i — 62) ■■■ 

x Sl:: 3 J::* (0„_j - 0n-2)A < 2 ■ ■ ■ A^A^. (B29) 

For this to happen the following consistency condition 
must hold 

n-2 

i>j 

n—3 

X ( _fj Pli^iP^ i li-\)P^n-\tn-Q.P(- l n _ 2 1n-3 

i=2 
n-2 

X J] Pe^Pe^ (B30) 
i=l 

This consistency condition is equivalent to what we need 
to establish in Ecrn. (|B26|) provided 

n-2 n-2 

II *Wi = II P <*n-i (B31) 

i=l i=l 

If the {ej}^ are all kinks or anti-kinks then {e^}^ is 
but a permutation and the above equality certain holds. 
If there are kinks and anti-kinks, then a kink- anti-kink 
pair a — a in {ejj^Zf may become /? — /3 in {e-}"^ 2 but 
because 

i^-i^Wi = Pffe^P^, (B32) 

the above equality is unaffected. 

5. Periodicity Axiom for Semi-local Operators of 
Section IV F 

We arc now able to write down the correct form for the 
periodicity axiom for the two-particle form-factors of the 
Sections III B 4. and 5. The periodicity axiom for the 
form factor of the scalar operator transforming as the 10 
(Section III B 4.) in the C-symmetric basis is 

f[abc] A {Q _ 2m) = f[abc ]A { _ 0) (333) 

The semi-locality phase in this case is 1. As this operator 
couples to two identical kinks, transforming to the non-C 
symmetric basis leaves the equality unchanged: 

f^ cU (0 - 2vri) = ftt c]A (-0) (B34) 
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We now consider the periodicity axiom for the form 
factor of the scalar operator transforming as the 6. For 
specificity we consider the operator 



01 



,§(-02 + 03+04) 



(B35) 



This operator couples to two same chirality kinks. In the 
C-symmetric basis that periodicity axiom reads 



2m) 



-ff3a(~ ) for a, /3 even 
for a, (3 odd 



(B36) 



We see that for this particular operator the even kinks 
coupling to the field have a non-trivial OPE whereas the 
odd kinks coupling do not. Transforming to the non-C 
symmetric basis introduces a sign changes as P a pPpa = 
-1: 



for a, P even 
for a, (3 odd 



(B37) 



If we had instead considered an operator of the form, 
ipi = e § (-02-03+04) the role of even and odd kinks would 



be swapped. But regardless of which particular scalar 
operator transforming as the 6 we consider, one type of 
kink will carry an additional sign in the periodicity axiom 
relative to the other. 



APPENDIX C: DERIVATION OF THE 
SUSCEPTIBILITY OF NON-INTERACTING 
LADDERS 

The non-interacting susceptibility is given by 

X° n (qx,uj) = X22{<lx,^) 

3 

= --(Li(u>,q x ) + L 2 (uj,q x ) + L 3 (u, q x ) + L±{lu , q x )) ; 
3 

= --(Li(u,q x ) + L 2 (uj,q x ) - L 3 (uj,q x ) - Li{u,q x )), 



(CI) 



where the L, 's are defined as follows: 



L 2 (u>,q x ) 



1 1 / 1 Atsm(q x /2)cos(K Fb + q x /2) , 4tsin(g x /2) cos(K Fb - q x /2)\ 

; tann tann ; - : 

2n T{uj, q x ) \ r{u,q x ) F{u,q x ) J 

11/ tanh -1 4t sin (g^/ 2 ) cos{K Fab + q x /2) _ , , -i 4i sin(g x /2) cos(K Fab - q x /2) 



iuj-\-e 



2nT(u,q x )\ T(u,q x ) Y{u>,q x ) 

L 3 (^q x ) = J dT&e^-^T^ac^ 

= L 31 (u,q x ) + L 32 (uj,q x ) + L 33 (u,q x ) + L 34 (u,q x ); 

if, aM . , , , , s -i2icos(gi) + z2tcos(<?i -q) 

L 31 (u!,q x ) = — / dq 1 0(2tcos(q 1 )+t±) , — — — j 

47T J (—iw + e — i2t±) z + 4i z (cos(<7i) — cos(gi + q)) £ 

= ^L 1 (u) + 2tj_,q x ) 

j i \ i [ a aio+ i \ + ^ -i2icos(gi) +i2tcos(q 1 - q) 

L 32 (u>,q x ) = — / dq 1 0(2tcos{q 1 )-t ± )-—— — — ■; — — 

Air J [—iw + e — i2t±) z + 4i z (cos(<7i) — cos(gi + q)) £ 

1 , 



2 

L 33 {oj 



1 / 

-L 2 (lo + 2tj_,q x ) 



\ i f /i/r, r \ \ iw + i2t± 

,q x ) = - J d qi 0(2tco S ( qi ) + t± ) ^ - - _ at±)2 - 4t2{cQs{qi) _ cQs{qi - q))2 

i 1 f _ lf T(LJ + 2t ± ,q x )t!w(K Fb + \q x \/2' 
-^T{u + 2t^q x )Y nK ^ /2 + i&n ( Hw _ i2 tx + c) 1 
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_! ^r(w + 2U,qx)tan(|K Fb - |q x |/2| 

+sgn(K Fb - |qx|/2) Trni i + tan 



|K Pb -|<u|/2| v - i2t ± + e) 



L 3i (u,q x ) = y dqi0(2tcos(qi) - tj_) 



(—iw + e - i2ijj 2 + 4t 2 (cos(gi) - cos(<?i + <?)) 2 



1 r _j r(cj + 2f_L, fe )tan(A> a6 + |g x |/2), 

""^Fab+lfol/a+ta 11 (" 



, ^ \ -i / r(co + 2t ± ,q x )tan(|K Fab - |q x |/2|) , 
+sgn(K Fab - |q x |/2)^n |KFab _ |qx|/2| +tan ( _ i2t± + e) . 



£4(w,&) = J dr^e i ^-^ K (T^(^a4 CT (a ;) T) CafeCT (a ;! r))(^ CT 4(0,0)c 6CT (0 ! 0)) 



= L 3 (Lu-tt ± ,q Xl K Fb <-> K Fab ), (C2) 

I 

where T(uj,q x ) = (-(w + ie) 2 + 16i 2 sin 2 ^^)) 1 / 2 and 
n g is an integer minimizing |?i g 7T — q\. 
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